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Motivated by the motion of biopolymers and membranes in solution, this article presents a for- 
mulation of the equations of motion for curves and surfaces in a viscous fluid. We focus on 
geometrical aspects and simple variational methods for calculating internal stresses and forces, 
and we derive the full nonlinear equations of motion. In the case of membranes, we pay particular 
attention to the formulation of the equations of hydrodynamics on a curved, deforming surface. 
The formalism is illustrated by two simple case studies: (1) the twirling instability of straight 
elastic rod rotating in a viscous fluid, and (2) the pearling and buckling instabilities of a tubular 
liposome or polymersome. 
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I. INTRODUCTION 

By definition, soft matter is easy to deform. Some ma- 
terials are soft because they have a high degree of symme- 
try. For example, the absence of long-range translational 
order in a liquid means that there is no resistance to 
shear in static equilibrium. Other materials, such as col- 
loidal crystals, have low symmetry but are soft since the 
fundamental constituents are large and the interactions 
between these constituents are of entropic origin, and 
therefore weak. For example, to estimate the shear mod- 
ulus 67 of a colloidal crystal, we assume a characteristic 
length scale a « 1 /xm and a characteristic energy scale of 
about ten times the thermal energy at room temperature, 
and find G w 10fc B T/a 3 « 4-1CT 2 Pa. This estimate is far 
smaller than the typical elastic modulus for a molecular 
crystal, where covalent bonds and molecular sizes deter- 
mine the size of the shear modulus. Finally, many objects 
are soft since they are thin in one or two dimensions. It 
is a matter of common experience that bending a thin 
rod or plate is much easier than stretching. Examples of 
thin objects at the micron scale abound: bacterial flag- 
ella, actin filaments, and cell membranes readily deform 
in response to thermal fluctuations or viscous forces aris- 
ing in fluid flow. The key to understanding the physics of 
these systems often hinges on understanding the dynam- 
ics of their shapes. Thus, the student of soft matter is led 
directly to the problem of describing how curves and sur- 
faces deform under external forces. Differential geometry 
is the natural tool for this task. In this review, we give 
a self-contained introduction to differential geometry of 
curves and surfaces that evolve in time. Our intended 
audience is someone who wants to use these equations to 
elucidate physical phenomena. Therefore, although our 
main focus is on mathematics, we give physical interpre- 
tations of the mathematics whenever possible. 

It is natural to ask why we need another review of dif- 
ferential geometry. There are introductory books of a 
purely mathematic al nature that are excellent starting 
point s for stude nts dCoxeterl. Jl969t|Millman and Parkeri 
119771 : iMorganl 119981: IStruikL 119881) . There are useful 
books with a broader mathematical scope, such as that 
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FIG. 1 A two-dimensional being [after iBurger and Starbirdl 
(2005)] considers the relation between the circumference and 
radius of a circle in his curved world; A" is a number he can 
determine by making careful measurements of length. 



of ( Nakaharal 120031 ). Finally, there are books and re- 
views that introduce diff erential geometry in th e con- 



text of solid mechanics (| Green and Zernl . 1968), fluid 



mechanics dAril. I1989D. and soft matter ( Davidl . 1989; 
Ivan Hemmen and Leiboldl . [20071 : IK amicn, 2002) . Our re- 
view is closest in spirit to these last four, but it is com- 
plementary since it emphasizes the variational point of 
view for deriving forces due to internal stresses acting on 
a stiff polymer or membrane. Like the review ( Kamienl 
2002). we give a unified treatment of curves and surfaces 
from the "extrinsic" point of view, in which distances 
along curves and surfaces sitting in three-dimensional 
space are measured using the notion of distance in that 
three-dimensional space. The alternative "intrinsic" ap- 
proach to the differential geometry of surfaces is to work 
only with quantities that can be measured by a two- 
dimensional being living on the surface. For example, 
the curvature K of a surface can be defined by compar- 
ing the circumference C of a circle on the surface to the 
radius R of the circle, as measured by the inhabitant of 
the surface. (See Fig. [TJ the formula relating K, C, and 
R displayed in the figure can be readily confirmed for 
the case of a sphere). Since polymers and membranes 
interact and bend in three-dimensional space, it is im- 
possible to formulate their equations of motions in terms 
of intrinsic quantities only, and we will not hesitate to 
characterize curves and surfaces with extrinsic quantities 
such as length in three dimensions. 

In addition to the tension between the extrinsic and 
intrinsic points of view, the issue of choosing coordi- 
nates always arises when studying curves and surfaces. 
Since geometric quantities such as vector fields and cur- 
vatures exist independently of the choice of coordinates, 
it is natural to defi ne these quan t ities w ithout invoking 
coordinates at all ( Misner et all fl973i) . Although this 
approach is elegant, it requires too much mathematical 



machinery for our purposes. In the case of curves, the ba- 
sic geometrical ideas can be captured with vector calcu- 
lus. Surfaces are more complicated and therefore require 
us to introduce some of the elements of the calculus of 
tensors, which we attempt to do with a minimal amount 
of formalism. 

A second reason not to insist on a coordinate-free ap- 
proach is that physical considerations often single out a 
coordinate system or a class of coordinate systems as spe- 
cial. For example, we will see that arclength parametriza- 
tion is natural for inextensible rods, and that convected 
coordinates — coordinates that are carried along by mate- 
rial particles — are natural for fluid membranes. A second 
point related to choice of coordinates is the distinction 
between an interface between two different phases and 
the physical surfaces that we study . For example, since 
a solidification front ( Langerl . Il980f) advances through a 
substance while the material points of the substance re- 
main at rest, we are free to use choose coordinates such 
that the in terface velocity is a lways normal to the tan- 
gent plane ( Brower et a/.l . 1 198*41 ). In contrast, the velocity 
of a material point in a polymer or membrane may have 
a component along the tangential direction. For motion 
of a polymer or membrane arising from deformation, this 
velocity component has physical consequences and can- 
not be removed by a change of coordinates. 



II. CHARACTERISTIC LENGTH AND TIME SCALES 
A. Thermal fluctuations and the persistence length 

To keep the review of manageable size, we limit the 
scope by disregarding thermal fluctuations. This as- 
sumption is a vast simplification, since the flexibility 
of polymers and membranes ensures that thermal un- 
dulations are always present. However, there are many 
shape-evolution problems in which these effects are small. 
Furthermore, the formalism described in this review 
is also useful even for problems where thermal effects 
play an important role, such as the dynamics o f ten- 
sion in semiflexible polymers ( Hallatschek et all , 120071) 
and t he surface tension of an und ulating fluid mem- 
brane ( Fournier and Barbettal l2008r i. In this section we 
briefly delineate the conditions under which thermal fluc- 
tuations are important. 

First consider a filament in solution. The filament 
could be a semiflexible polymer such as DNA or actin, or 
the flagellar filament of a bacterium (|BoaJ, 120021 ) . These 
filaments resist bending, and may be described by an en- 
ergy 



A f 1 

E =2 R* ds > 



(1) 



where A is the bending modulus, with dimensions of en- 
ergy times length, and s is distance measured along the 
filament, and R is the local radius of curvature of the 
filament. We will describe the curvature of a filament 
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in much greater detail in i jlll.AI for now the reader can 
simply imagine that the thermal motion of the molecules 
of the solution causes the filament to flex back and forth 
between arcs of constant curvature. If we imagine length- 
ening the filament, then our experience with macroscopic 
rods tells us that the filament is easier to bend. That is, 
for fixed filament stiffness A, the same force applied per- 
pendicular to the end of the filament leads to a greater 
deflection as the filament gets longer. Thus, the longer 
the filament, the more it will bend spontaneously due 
to thermal fluctuations. We can estimate the length £ p 
of a filament at which thermal fluctuations become im- 
portant by equating thermal energy k B T with the bend- 
ing energy of a filament bent into an arc of one radian: 
k B T ~ A/£ p , or £ p ~ k B T/A. The length £ p is called 
the persistence length. At room temperature, the persis- 
tence length of DNA is about 50 nm; £ p is of the order of 
microns for actin, and milli meters for m icrotubules and 
bacterial flagellar filaments (jBoall . 120021 ) . The effects of 
thermal fluctuations are small when the length L of a 
filament is small compared to t p , L <C £ p . 

Unlike interfaces which are simply the boundary be- 
tween two phases, membranes are comprised of molecules 
that resist bending. Therefore, the shapes of membranes 
are also governed by a bending energy. As we review in 
j jlV.Al there are two radii of curvature at any point of a 
surface. Thus, the bending energy of a fluid membrane 
is 



2 J \Ri R 2 



dS, 



(2) 



where K is the bending modulus or stiffness, with dimen- 
sions of energy, and dS* is the element of of surface area. 
Just as for filaments, we expect that as a membrane gets 
larger and larger, the thermal motion of the solvent is 
more and more effective at bending the membrane. How- 
ever, since the units of k do not contain a length, it is 
not immediately obvious how to estimate the persistence 
length using dimensional analysis. Therefore, we must 
use a more detailed approach. 

Consider a membrane that is weakly perturbed by 
thermal motion and spanning a frame with dimensions 
L x L. Since the membrane is almost flat, it can be 
described as the graph of a height function h over the 
xy-plane. In i jiV.AI we will see that the energy in Eq. © 
may be approximated as 



E 



(V 2 hfdxdy. 



(3) 



Our strategy for quantifying the effect of thermal fluctua- 
tions will be to estimate the scale L at which fluctuations 
in the normal to the membrane, n sa z — V/i, become of 
the order of unity, or in other words, large enough to spoil 
the assumption that the membrane is almost flat. First, 
use the equipartition theorem to assign to each mode of 
the quadratic energy function (|3|) the contribution k B T/2 
toward the average energy. The simplest way to apply 



the equipartition theorem is to assume periodic bound- 
ary conditions and write the height as a Fourier series, 
h(x,y) = exp[iq- (x, y))hq/L 2 , where the sum is over 
all wavevectors q = (2irm/L, 2im/ L) with the integers m 
and n less than a maximum value m ma x ■ The wavelength 
corresponding to the maximum is the smallest length at 
which the continuum theory applies, and is a molecu- 
lar scale £ mo i such as the thickness of the membrane. 
Thus, the mean-square amplitude of the fluctuations of 
the height are 



k B T 
S 4 ' 



(4) 



Longer wavelength ripples are easier to excite. We can 
now calculate the fluctuations in the normal direction n. 
Since the membrane is almost flat, it is convenient to 
calculate the mean- square fluctuations of the dev iation 
Sn = n — z w — V/i ( Chaikin and Lubenskvl . 0.995): 



(Sh ■ Sh) 



E 

qq' 



q q 



L 2 ( ft qV 



d 2 q k B T 
(2tt) 2 nq 2 ' 



(5) 
(6) 



The integral in Eq. © runs over wavevectors ranging 
from the short-distance cutoff g m i„ ss 2ir/£ mo \ to the long- 
distance cutoff given by the scale of the frame, q m&K ~ 
2-k/L. Assuming for simplicity a disc-shaped domain of 
integration, we find (Sn ■ Sn) w k B T /(2ttk) log (L/£ mo \). 
The scale at which the thermal fluctuations lead to sig- 
nificant deformation of the membrane is the length scale 
L = £ p for which (Sn ■ Sn) as 1, or 



oiexp 



2ttk 
kWT 



(>) 



For typical values £ mo \ ~ lnm and k ~ l0k B T, the per- 
sistence length is extremely large. Thus, thermal fluc- 
tuations are often disregarded in studies of the shape 
dynamics of vesicles. 



B. Time scales for stretching, bending, and twisting of 
filaments 

To understand the physical mechanisms that govern 
the dynamics of filaments, we describe the characteris- 
tic time scales for stretching, bending, and twisting. To 
get the orders of magnitude of these time scales, it is 
sufficient to use scaling arguments and work in the limit 
of small deflection from a straight filament. The precise 
nonlinear equations of motion require the development of 
the geometric tools presented in the later sections of this 
review. 

We begin by recalli ng some basic defin i tions from 
the mechanics of rods (lLandau and Lifshitd . Il986l ) and 
viscous fluids ( Landau and Lifshitd . Il987[ ). Consider a 
straight rod of radius a, with long axis aligned along 
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z and under longitudinal deformation. The stress a zz , 
or force per unit area acting through a cross-section at 
z, is proportional to the strain for small deformations: 
<j zz = Yd z u, where u is the displacement of the cross- 
section of the rod originally at z, and Y is the Young's 
modulus with units of energy per volume. The net elas- 
tic force per length acting on the element of size dz is 
therefore given by na 2 Yd 2 u dz. Now suppose the rod is 
allowed to relax, with the only resistance to the relax- 
ation being the viscous drag of the solvent. 

To estimate the viscous resistance, we first argue that 
at the small scales of polymer and cells, the typical vis- 
cous forces are much larger than the forces that are re- 
quired to change the velocity of the fluid. Viscous stresses 
are given by the product of the viscosity rj and the strain 
rate. Thus, the viscous force on an object of size £ is 
approximately 



iv 



( V v/£)£ 2 = r\vl, 



(8) 



where rj is the viscosity and v is a characteristic speed 
of the flow. The viscosity of water is 77 ~ 10 _3 N-s/m 2 . 
The product of mass and acceleration for an element of 
size £ of a fluid of density p is approximately 



Fi « (p£ 3 ){v 2 /£) = pv 2 £ 2 . 



(9) 



For micron-sized objects in water moving at speeds of the 
order of microns per second, the inertial force i 7 ! is about 
a million times smaller than the viscous force F v . This 
regime is the regime of "low Reynolds number," in which 
inert i a is irrelevant (lHappel and Brennerj . Il965at iHinchl . 
119881 : iKim and KarrilaL Il99lh . Thus, we may estimate 
the viscous drag on the element of the rod of length dz to 
be of the order of rjvdz. (In ^IIII we address the question 
of how the length scale a enters the drag.) 

We now have the ingredients required to estimate the 
relaxation time for longitudinal extension and compres- 
sion of a filament. Assuming that the inertia of the rod as 
well as the fluid is negligible, we balance the viscous drag 
with the elastic force, identify v with d t u, and drop nu- 
merical factors such as ir to find rjdtu ~ a 2 Yd 2 u. Thus, 
the characteristic relaxation time for longitudinal stretch- 
ing or compression of a rod of length L in a viscous fluid 
is 



l (L 

Y\a 



(10) 



To estimate the relaxation time for bending deforma- 
tions, note that the bending modulus A for a thin fil- 
ament of radi us a is related to the You ng's modulus 
by A ~ 7a 4 ( Landau and Lifshitzl . Il986f) . Just as we 
did for membranes, we can describe small deflections 
of a filament by a height function h(z), and approxi- 
mate the elastic energy as E w {A/2) J(d 2 h) 2 dz. By 
the principle of virtual work, the elastic force per unit 
length due to bending is given by the functional deriva- 
tive —5E/8h ~ —Ad z h. Once again disregarding numer- 
ical factors, we approximate the resistance to transverse 



motion as rjdthdz, and balance the elastic and viscous 
forces to obtain the bending relaxation time 



Tb 
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rj ( L 
Y 



(11) 



The ratio of the stretch relaxation time to the bend re- 
laxation time involves two powers of the aspect ratio, 
7s/ r b ~ ( a /L) 2 - This aspect ratio is typically very small, 
ranging from 1/50 for a 50 nm seg ment of DN A to 10~ 3 
for a typical bacterial flagellum (jBoall . l2002f ). There- 
fore, modulations in the extension relax much faster than 
bending deformations for polymer filaments. When the 
focus of interest is bending dynamics, thin filaments are 
typically assumed to be inextensiblc 

To estimate the relaxation time for twisting deforma- 
tions, suppose the filament is straight but twisted. We 
will precisely define twist for a filament with arbitrary 
contour in ^III.Al but for a straight filament define the 
angle 9(z) through which the cross-section at z has ro- 
tated due to twist. The torque applied through a cross- 
section is Cd z 9, and the net elastic torque on an element 
of length dz due to internal stresses is Cd z 9dz. 

The viscous torque resisting rotation at rate d t 9 is the 
product of the viscous stress r]d t 9, the surface area of the 
element ~ adz, and the moment arm a. Thus, the twist 
degree of freedom obeys the diffusion equation rja 2 dt9 = 



Cd z 9, and the twist relaxation time is 



r\a 2 l? 
C 



(12) 



For many materials, C ~ A ( Landau and Lifshitzl . 1 198"^) . 
and thus 



7] ( L 



T t ~ — - 

Y \ a 



(13) 



The twist relaxation time is much smaller than the bend 
relaxation time. Interestingly, the twist relaxation time is 
comparable to the stretch relaxation time. Nevertheless, 
it is common to consider the dynamics of inextensible 
elastic rods with both bending and twisting degrees of 
freedom, since twisting and stretching are usually driven 
by different physical processes. A simple example is pro- 
vided by an elastic filament, initially straight and twirling 
about its long axis ( £1111. F[) . Since the rotation induces 
twist but no stretching, it is convenient to treat the fila- 
ment as inextensible. 



C. Time scales for stretching and bending of membranes 

In this se ction w e revi ew the basic mechanics of 
membranes (ISeifertl . Il997ri . We confine our atten- 
tion to fluid membranes, in which there is no long- 
range order in the molecules making up the mem- 
brane. These membranes wrap up into closed sur- 
faces known as vesicles — liposomes in the case of lipid 
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molecu les, and polymer some s in the c ase of block copoly- 
mers ( Discher and Ahniedl . l2006t lYu and Eisenber j . 
Il998h . Both the lipids and the block copolymers are am- 
phiphilic, consisting of two parts with different affinities 
for water. The hydrophobic effect drives the molecules to 
form a bilayer, with the oily chains with the least affinity 
for water inside the bilayer. 

We now estimate the relaxation time for membrane 
stretching. Consistent with the discussion of ^II.Al we 
suppose the membrane is in the "high-tension" regime of 
stretching in which the resist ance arises from elasticit y 
rather than entropic effects ([Evans and Rawiczl Q990). 
Disregarding the bilayer structure for the moment, we 
may view the membrane as a two-dimensional fluid. Un- 
like the filaments we considered in £111. B| the membrane 
is a liquid, and to describe the resistance of the mem- 
brane to stretching and compression we use the density 
rather than a displacement variable. Let <p denote the 
number of molecules per unit area, and tf>Q denote the 
preferred number density. For small deformations, the 
clastic energy of stretching is 



E s = 



- 1 dS*, 



(14) 



where k is the compression modulus with units of energy. 
Since the modulus k is usually large enough that the den- 
sity is close to fa, work in terms of 5<f> — (f> — ef>Q. From di- 
mensional analysis, wc deduce that the two-dimensional 
pressure in the membrane, with units of force per length, 
is P2d = kS<p/4>Q. Gradients in the two-dimensional pres- 
sure lead to tangential forces in the membrane, which 
drive regions of non-uniform density to relax to constant 
density. The balance of these tangential forces with dis- 
sipative forces determines the stretch relaxation time. 
We consider the relaxation of a density modulation with 
wavenumber q. 

First we suppose that the dominant dissipative forces 
arise from the viscous forces the solvent exerts on the 
membrane. For a characteristic flow field v near the 
membrane, the stress is ~ r/Vw ~ rjqv. Thus the bal- 
ance of the gradient of pressure with the viscous stress 
leads to 



qkS(f>/4>o ~ m v - 



(15) 



Assuming the no-slip boundary condition between the 
solvent velocity and the lipid velocity holds, we may use 
V to represent the lipid velocity as well. The rate of 
change of the density 4> is related to the lipid flow by 
the continuity equation expressing conservation of parti- 
cles. In HTV.DI we will see how to formulate the continuity 
equation on a moving curved surface, but for our scaling 
argument it is sufficient to assert that dt<f> ~ V(0t> ), or 



54>/t s ~ q(j) v, 



(16) 



where r s is the stretching relaxation time of the mem- 
brane. Combining Eqs. (fT5|) and ([IB]) leads to 



rj_ 

kq ' 



or, taking the smallest wavenumber to be set by the size 
L of the vesicle, r s ~ rjL/k. 

To estimate the relaxation time for bending deforma- 
tions with wavenumber q, we balance the bending elastic 
force per unit area with viscous stress or traction from the 
solvent. We will give the precise form of the elastic force 
per area on a deformed membrane in £jIV.B.2[ but work- 
ing in analogy with our expression for the bending force 
per unit length on a filament, we have / ~ kV 4 /i ~ Kq 4 h. 
Again invoking the no-slip condition to identify dth with 
the velocity in the solvent near the membrane, we esti- 
mate the viscous traction on the membrane as ijqh/r^, 
where Tb is the bending relaxation time scale. Balanc- 
ing the_elastic_jOTce_rjer area with the traction leads 
to ( Brochard and Lennonl . [l975h 



Tb 



JL 

Kq 



3 • 



(18) 



Since the smallest q scales with vesicle size L, we may 
also write Tb ~ r/L /k. A simple scaling estimate for 
the magnitude of the bending stiffness is k ~ fc^ lol , 
where £ mo i is a molecular sca le such as the thickness of 
the membrane ( Seifertl . Il997f) . Thus, the ratio of the 
stretch relaxation time to the bend relaxation time for 
membranes scale as r s /rb ~ (lmoi/L) 2 , and just as in the 
case of filaments we expect density modulations to relax 
much faster than bending modulations. To estimate the 
time scales directly from the measured parameters for 
the lipid stearoyloleoyl phosphatidyl choline (SOPC), w e 
use k w 0.2 N/m and k m 1CT 20 J (jPimova et a/1 . l2002f) . 
which yields for a vesicle with L — 10 fim a stretch relax- 
ation time of t s w 10 -7 s and a bending relaxation time 
of r as 100 s. 1 Note that the bending relaxation time 
depends strongly on length scale; a ripple with a 1 fim 
wavelength relaxes 1000 times faster than a ripple with 
a 10 jitm wavelength. Thus, the time scale for bending 
dynamics of lipid membranes can easily be of the order 
of seconds. 

In addition to the viscous forces from the solvent, there 
are two other dissipative forces acting on the membrane: 
viscous forces due the surface viscosity of the membrane 
itself, and friction forces arising from the slipping of the 
two monolayers past each other. The relative importance 
of these forces changes with scale, with bilayer friction 
and then membrane vi scosity becoming dominant as the 
lengt h scale decreases (jSeifertl . Il997t ISeifert and Lange"rl 
1993). The implications of bilayer structure, both for 
dyna mics and for morphology, are described in ( Seifertl . 
I1997T) and will not be discussed in this review. How- 
ever, we will treat intrinsic membrane viscosity, since it 



(17) 



1 In principle, the relaxation dynamics of the compression modes 
could be so fast that inertial effects become relevant, invalidat- 
ing our estimates for viscous friction. This comment applies to 
both filaments and membranes, but does not change the conclu- 
sion th at the compression m odes relax much faster than bending 
modes l iSeifert and L anger. 199jj|). 
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FIG. 2 (Color online.) The directors for a rod are chosen 
along the principal axes of the cross-section of the rod. 



is even more important than solvent viscosity in the case 
of polymersomes. 

Surface viscosity relates the viscous force per unit 
length f v to the shear rate, / v ~ fiVti. Thus, the ratio of 
surface and bulk viscosity is a len gth, Isn = nj 1 7, term ed 
the "Saffman-Delbriick length" (jHenle et ali . l2008t l in 
honor of Saffman and Dclbriick's study of the mobility of 
particles embe dded in a thin liquid film atop a deep liq- 
uid subphasc (Saffman n and Delbruckl . Il975h . Since the 
force per unit area acting on an element of membrane is 
given by a derivative of the force per length, the surface 
viscous force per unit area scales as /iq 2 v compared to 
r\qv for the traction from the solvent. Therefore, surface 
viscosity effects become dominant when £sd9 > 1- Sur- 
face viscositi es for lipid membranes are of th e order of 
lCT 9 N-s/m (|Dimova et ali 120061 Il999l l2000h . whereas 
the surface viscosity of di block copolymer ves icles can 
be about 500 times higher (jDimova et al. I. I2002D . There- 
fore, the Saffman-Delbriick length is ^sd ~ 1 A*m for li- 
posomes, and £sd ~ 1 mm for polymersomes. Therefore, 
surface viscosity effects dominate the dynamics of poly- 
mersomes. However, we will see in CVTDl that solvent vis- 
cosity cannot be neglected in the dynamics of almost fiat 
membranes. Nevertheless, the bending relaxation time 
for a polymersome may be obtained by the same argu- 
ment that led to Eq. (fTS)) with fiq replacing rj, leading 
to 



Tb 



JL 

2 



(19) 



or Tb ~ [xL 2 /k for a vesicle of size L. Since surface vis- 
cous stresses will resist the relaxation of both compres- 
sion and bending modes, the relation r s <C Tb still holds 
for polymersomes. 



frame of directors. Then we use variational arguments 
to derive the elastic forces per unit length acting on an 
element of a rod in an arbitrary configuration. After a 
brief discussion of compatibility relations for rod strain 
and angular velocity, we conclude the section with an 
overview of the viscous forces acting on slender rods in 
the overdamped regime of Stokes flow. 



A. Rod kinematics 

The centerline of a rod is given by r(£), the position of 
the point with coordinate £. Since £ is fixed for a given 
material point, we call £ a material coordinate. Alterna- 
tively, £ is sometimes called a convected coordinate since 
the label £ is carried along with the material point as the 
body deforms. As mentioned earlier, it is much easier to 
bend a rod than it is to stretch it. Thus, we consider 
inextensible rods. Note that £ is a material or convected 
coordinate; it is a label for material points of the center- 
line of the rod. It is convenient to take the parameter £ 
to run from to 1 from one end of the rod to the other. 
The distance ds between two nearby points on the rod 
centerline is given by the Pythagorean theorem and the 
chain rule, 



(20) 



The vector c^r = dr/dt; is tangent to the centerline; 
normalizing this vector leads to the unit tangent vector 
r' = d s r. 

Our characterization of the configuration of a rod is 
completed by specifying the orientation of the rod cross- 
section at each s. Let di and d2 lie along the principal di- 
rections of the rod cross-section (see Fig. [2J . Then define 
d.3 = di x d2 . For simplicity, we will consider unshearable 
rods, for which the cross-section always remains perpen- 
dicular to the rod centerline, or 0I3 = r' ( Antman . ll995h . 
Note that if the cross-section is circular, then we must 
arbitrarily choose a direction di at each point s. For the 
case of a rod which is straight in the absence of external 
forces, a sensible choice is to make di(s) constant. 

Let the subscripts a, 6, ... = 1, 2, or 3 label the mem- 
bers of orthonormal frames. The configuration of the 
rod is given by the orientation of the frame for all s. In- 
stead of specifying the orientation of the frame, it is much 
more convenient to describe the configuration of the rod 
in terms of how fast the orthonormal frame rotates as s 
increases. The frames for any two nearby points differ by 
a small rotation. 



III. FILAMENTS 



d' = n x d„ 



(21) 



In this section we derive the equat ions of motion for 
flexib le filaments in a viscous liquid ( Doi and Edward"sl 
Il986l ) . We begin with kinematics and review how to com- 
pactly describe bending and twisting using a body-fixed 



where in the body-fixed frame of directors, fi = Qidi + 
fl2d.2+fl3d.3. We will use the Einstein convention of sum- 
ming repeating indices to save space: J7 = f2 a d Q . The 
infinitesimal rotations £l a are uniquely defined once the 
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FIG. 3 The basis vectors of the material and Serret-Frenet 
frames in the plane of the rod cross-section. 

directors are chosen, since f2i = d' 2 ■ d.3, and fl2 and O3 
are determined by cyclically permuting the indices. The 
quantities f2i and ^2 are the components of the curvature 
of the centerline, and is the rate at which the cross- 
section twists around d3 as s increases. For comparison, 
it is useful to recall the Serret-Frenet frame curve of a 
space curve, {T, N, B}, where T = r' is the unit tangent 
vector, and the normal N and the binormal B are defined 
by 

N' = tB-kT (22) 
B' = -tN (23) 
T' = kN. (24) 

The curvature k is the rate of rotation of T about B, 
and the torsion r is the rate of rotation of N about T. 
It is conventional to take k > to make N point to the 
center of curvature. Comparison of Eqs. (f2"Tj) and (|2"2Tl2"4"|) 
reveals that k 2 — O J -r-f2| an d ^3 = ip' + T , where tp is the 
angle between N and d.2: d.2 • N = cos?/' (Fig. Note 
that the directions of the normal and binormal become 
ambiguous when the rod is straight, k = 0. This ambi- 
guity can be sidestepped by using the so-called natural 
frame, which is well-defined ev e n at points where k = 
(|Bishod . Il975t iGoldstein efail , 1 19981: IWolgemuth et all 
120041) . However, there is no ambiguity about di and d2 
when the rod is straight; note further that fii and f^2 can 
have either sign. In this review, we consider rods which 
are straight in absence of external forces. Thus, the f2 a 
may be interpreted as strains associated with bending 
and twisting. Note that a rod with fl = K0CI2 is an arc 
of a circle, and a rod with tt = Kod.2 + Tgd^ is a helix. 
Note also that all three variables £l a must be given to 
specify the shape of a rod; specifying only the curvature 
and twist does not determine the shape. For example, 

^circle = 4 cos(s/i?)di + — sin(s/i?)d 2 + 4 d 3 (25) 

XI Xl XI 

describes a rod bent into a circle with curvature 1/ R and 
twist 1/R, whereas 

S^hdix = -5CI2 + -^d 3 (26) 




FIG. 4 Two rods with the same curvature k and the same 
twist fi3. For each rod, the black line traces the path of the 
tip of di. Note that the fact that one rod is closed and one 
is open is immaterial for determining the twist. 



describes a helix that also has curvature 1/R and twist 
1/R (see Fig. gj). 



B. Energies, forces, and moments 

Our next task is to find the force ids exerted on an 
element of the rod of length ds by the rest of the rod. 
This force equals the force the element exerts on the sur- 
rounding medium. It is useful to first consider a simple 
but classic question: what is the force per unit length 
that a curve in the plane exerts on its surroundings when 
its energy function is E = 7 J ds? Suppose also that the 
curve encloses a fixed area. Thus our problem models a 
two-dimensional droplet of incompressible fluid. To find 
the force for an arbitrary configuration r(s) of the curve, 
we impose an external force per unit length equal and op- 
posite to f (s) to put the instantaneous configuration in 
mechanical equilibrium. Thus, f is given by the principle 
of virtual work, or the condition 6E tot = 0, where 

SE tot =5E + J f(s) ■ 5r(s)ds - pSA, (27) 

and p is the Lagrange multiplier enforcing the constraint 
of fixed enclosed area A. To calculate SE tot , we must 
find how the length ds changes under a variation of the 
curve (Fig [5]): 

Sds = yj[r' + (<5r)'] 2 ds - ds 

= r' • (<5r)'ds, (28) 

where the second line of (|28p is accurate to first order 
in <5r. Note that r' = d 3 is the unit tangent vector. 
Our formula correctly shows that the length element is 
invariant under rigid motions, such as the translations 
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r(si) + Sr(si) 




FIG. 5 In general, the variation of a curve stretches the curve. 



Sr = a and the rotations Sr = 8 x r, where a and 9 are 
constant vectors. Using Sds to calculate SE to t yields 



SE tot = 7 / r' • (<5r)'ds + / (f-pN) ■ Srds 

-jr" + f-pN] -Srds, (29) 



where we integrated by parts and used the fact that 
Sr(L) = Sr(0) for a closed curve. Since Sr is an arbi- 
trary variation, the force is given by 



f = (jK+p)N. 



(30) 



Note that the curve can only support normal forces per 
unit length. In equilibrium, f = 0, and we recover the 
Young-Laplace law, with the Lagrange multiplier p iden- 
tified as the difference in the pressure outside the curve 
and inside the curve. Since p is constant, n is constant, 
and the length-minimizing curve is a circle. 

We now appl y this approach to inext ensible rods that 
resist bending (jHallatschek et ~al\. 120071) and twist. Un- 
like the previous example of a one-dimensional interface 
in the plane, an elastic rod has internal structure with 
its state defined by the strains f2 Q . Suppose the elastic 
energy of the rod is given by 



E = J £(fii,fi 2l fl3)ds. 



For example, the energy density 



21 2 3 2 2 3 



(31) 



(32) 



describes a rod with bend modulus A, twist modu- 
lus C, and a straight shape in the unstressed con- 
figuration. Our formalism is easily generalized to 
handle m odels in which £ de p ends on derivatives of 
the n a (|Kan and Wolgemuthl . |2007t) . such as mod- 
els for pol ymorphic phase transformations in bacte- 
rial flagella (ICoombs et all l2002HGoldstein et all [2005 
ISrigiriraiu and PowersL I2005L l2006t) . However, we will 
disregard such "strain-gradient" terms for simplicity. 

It is natural to suppose that the ends of the rod are 
subject to external forces F cxt (0) and F ext (L), and mo- 
ments M cx t(0) and M cxt (L). To hold an arbitrary con- 
figuration in equilibrium with these boundary conditions, 
we must apply an external force per unit length — f(s). 
Since the rod can twist about the tangent d3, we must 



d 3 + M 3 



do 



d 2 + Sd 2 



FIG. 6 Rotation of directors under a variation with S\i 7^ 0, 
<5X2 = Sxi = 0. 



also apply an external moment per unit length — 77i(s)d3. 
As in the closed curve example, f (s) and m(s)A^ are the 
forces and torques per unit length that the rod exerts on 
the medium for a given configuration. These are again 
given by the condition SEtot = 0, where SE io t includes 
the work done by all external forces. Note that we are not 
considering the most general form of the external torque 
per unit length. A torque per unit length parallel to the 
centerline of the rod arises naturally when considering 
rotation of a filament in a viscous liquid (see COLF) . but 
our formalism can be easily modified to handle exter- 
nal torques per unit length that are perpendicular to the 
centerline. 

Since the rod is inextensible, care must be taken when 
performing the variation, since most variations Sr(s) 
stretch the rod (Fig. [5]). The three components of 
Sr are not independent for a variation that conserves 
the distance between material points on the rod. It 
is convenient to introduce a Lagrange multiplier func- 
tion enforcing r' • r' = 1 and then vary the energy 
with r espect to the three compon ents of Sr indepen- 
dently ( Goldstein and Langerl . [l~995( ) . Instead we will fol- 
low the less common procedure of dispensing with the 
Lagra nge multiplier func t ion a nd constraining the vari- 
ation (Hallatschek et al. 

1, 12007ft . Since the constrained 
variation does not stretch the rod, arclength is conserved, 



Sds = 0. (33) 

Therefore, the order of applying the variation 5 and 
derivatives with respect to s is immaterial: 

5d s = d s S. (34) 

We may drop the distinction between S(r') and (Sr)', and 
simply write Sr' . 

Since the energy depends on the tt a only, the configu- 
ration of the rod is completely determined by the orien- 
tation of the directors in sp£iC6, d a . Thus, we may vary 
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the configuration without stretching the rod by subject- 
ing the orthonormal frame of directors at each point s to 
an infinitesimal rotation, <5d a = e x d a where e = <5x Q d a . 
The infinitesimal angles 5x a are related to the directors 
via 



Since the energy is a function of f2 a only, we have 



Sd? ■ d.s 



(35) 



with 6x2 an d 5x3 determined by cyclically permuting the 
indices in (|35|) . Thus, under a variation, the orthonormal 
frame vectors d2 and d3 rotate about di through the an- 
gle Sxi (Fig. O. The two angles S\i and 8x2 correspond 
to the two degrees of freedom of Sr under the constraint 



&Xi = -d 2 ■ (St)' 
5 X 2 = di • (5r)'. 



(36) 
(37) 



In general, the are no globally defined angles Xa whose 
variations correspond to the infinitesimal angles Sxa- 
The exceptions are cases when one of the directors is 
constant for the whole rod. For example, we may take 
di to be normal to the plane containing a planar curve; 
then Xi ma y be defined as the angle between the local 
tangent vector and a fixed direction in the plane, and x'i 
is the curvature. Likewise, the director d3 has a constant 
direction for a straight, twisted rod, with X3 the cumu- 
lative angle through which di has rotated, and x'3 the 
twist Q3. 

The variation in the strain is found by varying Eq. (|21[) , 
using Eq. (CM)) , 5d a — ex d a , and e = 5x a d a , and finally 
noting that V x d a = for all a implies V = 0: 



sn 



e x fi. 



(38) 



It is also useful to write the variation in component form: 

sn, = S X [ + n 2 5 X 3 - n 3 5 X 2 (39) 
SQ 2 = S X ' 2 +n 3 S Xl -Q 1 S X3 (40) 
5Q 3 = 6x3 + ^5x2-^25x1- (41) 

The components of O x £ and not e x $1 appear on the 
right-hand sides of Eqs. ([59THT]) . since the basis vectors 
d a in 5Q and e must be varied and differentiated, respec- 
tively. To interpret these formulas, note that a variation 
with e = <5x3(s)d3 changes the twist Q3 and the princi- 
pal curvatures fii and ^2, but not k = y/Clf + f2|. A 
deformation of the centerline, e = <5xi(s)di + 5x2(s)d2 
changes the twist if the rod is curve d. This fact under- 
lies the coupling of twist and writhe rtCalugareanul. 11959) : 
Dennis and HannavL 200 5: FullcrL ll971t iGoldstein et all . 
19981 : iKamienL 1 19981: IWhitd . Il969h . 

We now have all the necessary ingredients to find the 
force per unit length f . We demand 5E to t = 0, where 



5E U 



5E 



({■Sr- 



- F ext (L) ■ 5r(L) 

- M ext (L) ■ e(L) - 



1128x3) ds 

- F ext (0) • Sr(0) 
■Mext(0)-e(0). (42) 



SE 



M a 5Q a ds, 



(43) 



where M a = d£/dil a - The formulas for the variation of 
the strain, Eqs. (|39H41|) . imply that 



Mjn a = M a 5 X ' a + M • O x e (44) 
= (M a 6 Xa )' - M' a 8 Xa - e ■ n x M (45) 
= (M a S Xa y-M'-e. (46) 

Likewise, if we introduce the force F(s) acting on the 
cross-section of the rod at s by the relation F' = f , we 
have 



F'-St 

(F • Sr)' — F- St' 
(F • St)' - e • d 3 x F. 



f ■ Sr = F' ■ Sr (47) 

(48) 
(49) 

Combining Eqs. ([46]) . ([49]) . and (|42]) . and demanding 
5E to t = for arbitrary (inextensible) <5r yields F(L) = 
F ext (L), F(0) = -F ext (0), M(L) - M cxt (i), M(0) = 
-M cxt (0), and 



M' + ds x F. 



(50) 



The boundary conditions on M lead us to identify M(s) 
as the moment exerted through the cross section at s by 
the part of the rod with arclength greater than s on the 
part with arclength less than s. Note that when m = 0, 
Eq. (|50D i s the equilibrium momen t balance equation for 
a rod (jLandau and Lifshit d . Il986l ). 

For the rod energy density (|3l?l) , the relation M a — 
d£/dn a yields M = A(f2idi + fi 2 d 2 ) + Cn 3 d 3 , or 

M = Ad 3 x d 3 + Cn 3 d 3 . (51) 

The transverse part of Eq. ([50]) yields 

F = -At'" + Cn 3 r' x r" + Ar'. (52) 

Our variational procedure does not determine the func- 
tion A. The value of A is chosen to enforce the con- 
straint r' • r' = 1, but from Eq. ([52]) we can s e e that 
it may be interpreted as a ten sion [see iDeutschl ( 19881 ) 
or IGoldstein and Langerl ( 19951 ) , for example] . Note that 
—Ar'" has a tangential part. 

To summarize this section, we have used the principle 
of virtual work to derive the moment M and force F act- 
ing through any given cross-section of a rod with energy 
density ([32]) . The constitutive relation M a = d£/dft a 
arises naturally from the variation. Likwise, the condi- 
tions of force balance F' = and moment balance M' + 
d3 x F ~ in equilibrium, when f = and m — 0, are not 
imposed but are implied by the principle of virtual work. 
These statements can be generalized to arbitrary energy 
densities a long a curve (jStarostin and van der Heiidenl . 
I2007L I2009T) . 



10 



C. Compatibility relations for rods 

The variational formula Eq. (f38|) leads to important 
relations between the angular velocity and strain. The 
local angular velocity vector uj is given by w = e/5t as 
St — ► 0, where the variation is interpreted as the change of 
the rod shape over a time interval St. Thus Sd a = e x d a 
implies 



ddg 

at 



■■ ijj x d„ 



(53) 



Likewise, since we may similarly define dfl/dt = 6Sl/St 
asSt^O, Eq. ([35]) implies 



~dt 



(54) 



These compatibility relations arc also conveniently de- 
rived using th e rotation matrices tha t act on the di- 
rector frames ([Wolgemuth et al\ , l2004h . A third way of 
deriving the compatibility relations is to recognize that 
since s and t are independent variables, the order of tak- 
ing derivatives does not matter, and thus we must have 
d 2 d a /dtds — d 2 d a /dsdt. We will see below that the 
3-component of Eq. (|5T 



dt 



du> 3 
ds 



(55) 



is required to obtain the evolution equation for twist. 
Once again, it is useful to note that uj' contains deriva- 
tives of iv a and derivatives of the directors when compar- 
ing Eq. and Eq. ([M]) ; compare with the discussion 
after Eqs. (ptyllj) . 

We may use Eqs. (|2Tj) and (|53|) to recast Eq. ([55]) in 
terms of the shape of the filament r: 



9A = ^ 3 + r'xr".(9 t r)'. 



(56) 



This relation may be viewed as a conservation law 
for twist density ^3, with L03 the twist current, and 
the second term on the right-hand side acting as a 
source for twist (I Goldstein et al\ . \l99& iKamienl . [l99l 
iKlapper and Tabon 1994h . The first term may be un- 
derstood by observing that the only way the twist in a 
straight rod can change is if the rate of rotation W3 is 
not uniform along the rod. The second term reflects the 
fact that motions of the backbone that do not involve a 
rotational velocity can also change the twist. Such 
motions require a change in the writhe, where the writhe 
of a closed curve is computed by counting the number of 
self-crossing of the curve in a planar projection, and then 
averaging over all projections [see iDennis and Hannavl 



(2005) and references therein]. Although writhe is de- 
fined for a closed curve, its variation has a local form 
and can be considered for an open curve dGoldstein et ali , 
119981: iKamienl . Il998t IKlapper and Tabori . Il994h . 



D. Viscous forces: slender body theory and resistive force 
theory 

In this section we describe how to approximately cal- 
culate the forces acting on a slender body immersed in a 
viscous fluid, and give a slightly more quantitative discus- 
sion of viscous fluid dynamics tha n in flll.Bl The reader 
is directed to ( Lauga and Powersl l2009h and references 
therein for a more detailed treatment. The stress in 
an incompressible viscous fluid is related to the strain 
rate by <7jj = —pSij + rj^iUj + Vjttj), where p is the 
pres sure, 77 is the visco s ity, a nd u is the flow veloc- 
ity (jLandau and Lifshitzl . 1 19871) . Note that we use the 
indices i, j, ... for the Cartesian coordinates. As argued 
in ijll.Bl viscous forces dominate at the small scales of 
polymer and cells. Thus we are permitted to disregard 
inertia. In this limit, the flow is governed by the Stokes 
equations: 



?7V 2 u - Vp 
Vu 





0. 



(57) 
(58) 



Note that the pressure p plays the same role as A did 
in our discussion of inextensible polymers, Eq. (f52"| : the 
value of p is determined to enforce the constraint of in- 
compressibility, V • u = 0. The boundary conditions for 
the Stokes equations are that the velocity of the fluid at a 
solid boundary must match the velocity of the boundary. 

Ideally, we would like to calculate the viscous forces 
acting on a moving and deforming filament. Since no 
time derivatives appear in the Stokes equations (|57H58p . 
the flow u and thus the forces acting on the filament 
at a given instant depend only on the geometry of the 
filament — its velocity and shape at that instant. How- 
ever, it is generally impossible to find analytic expres- 
sions for the forces for anything but the simplest shapes. 
The problem can be simplified by developing an approx- 
imation scheme that exploits the smallness of the ratio 
of the rod diameter to its length; this ap proach is known 
as slender-body theory ( Lighthill Il976l ). Slender-body 
theory relies on the linearity of the Stokes equations to 
construct solutions for the flow around a moving or de- 
forming filament by superposition of singular solutions. 
These singular solutions are the flows induced by a point 
force or a dipole source, and are analogous to the so- 
lutions used in the method of images in electrostatics. 
The end result is that slender-body theory provides a 
method, typically numerical, for computing the hydro- 
dynamic forces distributed along a fil ament. Repr e senta - 
tive examp l es inc lude work by I Shelley and Uedal ( 20001) 
and lCorted (|200lh . 

In the limit of an asymptotically thin filament, slender- 
body theory simplifies to a local theory known as resistive 
force theory (|Chwang and Wul . ll9Tlb iGrav and Hancock! 
Il955l ; Lighthill. Il976l ). Resistive force theory is local in 
the sense that it describes the hydrodynamic force per 
unit length acting at a point on a filament in terms of 
the velocity of the filament at that point, in contrast with 
slender-body theory, which accounts for the fact that the 
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motion of every point along the filament induces a flow 
that affects the hydrodynamic forces at any given point 
on the filament. Thus, it is often stated that resistive 
force theory does not properly account for hydrodynamic 
interactions. Nevertheless, since resistive force theory is 
much simpler to implement that slender-body theory, and 
because it can give accurate results, we will use resistive 
force theory to formulate the equations of motion of fila- 
ments. 

In resistive force theory, the fluid is treated as a passive 
background material that does not respond to the motion 
of the slender object. The viscous force per unit length 
acting on the rod is anisotropic, and given by 



f v = -Cl[v - ( d 3 • v)d 3 ] - C||(d 3 • v)d 3 



(59) 



where v is the velocity of the rod relative to the local 
velocity of the fluid. In the limit of small a/A, the friction 
coefficients in Eq. (15^1) are defined by 



Cii 



47r^/ln(A/a) 
27r/i/ln(A/a), 



(60) 
(61) 



where a is the filament diameter, and A is a cutoff rep- 
resenting either the rod len gth or the characteristic ra - 
dius of curvature of the rod ( Ke ller and Rubinowl . ll97'6h . 
Note that replacing A by 2 A in Eqs. ([501 EU merely leads 
to subdominant additive terms in the denominator. The 
expression (|59p becomes accurate when a becomes much 
smaller t han the rod's length or ty pical radius of cur- 
vature A ( Keller and Rubinowl . 1197a ). Macroscopic scale 
experiments show that resistive force theory does sur- 
prisingly well in describing the sh ape of driven elastic 
filam e nts in high l y visc o us fluid ( Koehler and Powersl 



l2000t iQian et all , 120081: \Yuet allBoO^ . Numerical 



calculations show that resistive force theory is inac- 
curate for tightly coiled helices in flow, where hydro- 
dyna mic interactions are ex pected to be more impor- 
tant ( Kim and Powers! . 12005). Even in a quiescent fluid, 
resistive force theory can give qualitatively wrong pre- 
dict ions for the sense of rotation of a towed superhe- 
lix (jjung et al\ . 120071 ). Thus, resistive force theory is 
usually justified in a first approach to a new problem, 
but its predictions must always be viewed with a critical 
eye. 

In some situations, suc h as the rotation of a n early 
straight rod about its axis (jWolgemuth et aUl2000| ). it is 
necessary to include the viscous moment per unit length 
for rotation about the local tangent vector d 3 . The ap- 
propriate local approximation for this moment is the vis- 
cous torque per unit length on a rotating cylinder: 



-Cr(d 3 



(62) 



where C 



Anria 2 is t he r otational drag coeffi- 



cient ( Landau and Lifshitd . Il987l ). In other situations, 
such as the rotation of helix, the contribution to the 
torque due to the rigid body motion of the centerline 
dominates the contributions from Eq (|62p . 



E. Equations of motion for curves 

We are ready to combine all the parts we have intro- 
duced and arrive at the coupled equations for the evolu- 
tion of a rod. Since we work in the highly overdamped 
regime, all forces must sum to zero: 



F' + f v = 0. 



(63) 



Suppose that there is no imposed flow in the fluid; thus 
the velocity of the rod relative to the fluid is v = dtY. 
Using the elastic force on a cross section (|52[) and the 
viscous force per unit length (|59p . we find 

a(d t r) x +C|| (d t v) n = -Ar"" + C(n 3 v' x r")' + (Ar')', 

where (d t r)± = d t r - [d 3 • (<9 t r)]d 3 and (d t r)\\ = [d 3 • 

(<9tr)]d 3 . To determine A, write v in terms of the mo- 
bilities £j_ and 0| and the internal force per unit length 



F', 



C|| Cx' 

and use d t (r' • r') = 2d t (r') • r' = 29 s (v) • r' to find 



d J ■ d 3 



(65) 



(66) 



Using Eq. ([52^1 to calculate f = F' yields the desired 
differential equation for A. 

The elastic and viscous moments must also sum to 
zero: 



M' + da x F + m v 



0. 



(67) 



The only nontrivial terms in Eq. (|67p are the tangential 
components, which simplify to 



cn' 3 = r> 3 . 



(68) 



To find the evolution equation for twist, combine 
Eqs. and flS| to find 



Crfltfia = Cn' 3 ' + C r r' x r" • (5 t r)' , 



(69) 



where dtY must be deduced from Eq. 

We close this section by recalling the characteristic 
time scales for the motion of filaments in a viscous 
fluid described in ijlLBl Examination of Eq. ([51)1 for a 
nearly straight filament with no twist reveals that bend- 
ing modes have a characteristic relaxation time Tb ~ 
(±L 4 /A. Likewise, Eq. (|6"9"]) applied to a straight but 
twisted rod implies that twist deformations have a char- 
acteristic relaxation time T t = ( T L 2 /C. 



F. Illustrative example: twirling and whirling of elastic 
filaments 

To illustrate the use of the equations of motion, we con- 
sider a fundamental problem: the whirling instability of 
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FIG. 7 An elastic filament immersed in a liquid with viscosity 
r), and twirled about the z-axis with angular speed too. The 
straight state is unstable above a critical angular speed. 



rod writhes into a nonplanar shape once it is twisted by 
about one turn. In our problem with distributed exter- 
nal moments, the speed at which the typical twist torque 
balances the characteristic bending torque is given by 
CrWcriti ~ A/L, or 



Wcrit ~ A/i^L 2 ). 



(71) 



Since C » A, the critical rotation rate is the same order 
of magnitude as the relaxation rate for twist. 

To find the precise value of the critical speed, linearize 
the equation of motion about the base state, Eq. 



a rotating elastic rod i n a viscous fluid dLim and Peskinl . 
12004 1 Wada and NetzL |2005 IWokemuth et a/.l . l2000h . A 
motor at s = rotates the rod about its long axis 
with angular speed loq (Fig. [7|). Initially, the rod is 
straight. Viscous forces cause twist to build up. For 
sufficiently high rotation speed, the straight state is un- 
stable, and the rod writhes. Study of this problem will 
provide us with intuition that will help us understand 
more complicated phenomena in rotating filaments, such 
as the propag ation of polymorphi c instab i lities in bac- 
terial flagella (jCoombs et al 1 120021: lHotanlll982f) or the 
chiral dynamics of beating cilia in embryonic develop - 
ment (|Hilfinger and Julicherl . 120081 : iNonaka et a/1 . 119981 ) . 

As mentioned earlier, in this review we disregard the 
effects of thermal fluctuations. We consider a rod with 
length much smaller than the bend persistence length 
and twist persistence length. We also suppose the length 
is much smaller than the twist persistence length k^T /C, 
which is defined by an argument similar t o the one we 
gave in ClA"l See (jWada and Netzl . 120061) for a study 
of the whirling instability for rods so long that thermal 
effects are important. 

Since we assume the rod is thin, it is natural to use the 
local drag formulas Eq. (|59")) and Eq. (|rj2")l to describe the 
viscous forces and moments acting on the rod. Dropping 
the subscript 3 to reduce clutter, the balance of torques 
about the tangential direction ([rI5)) becomes CO,' = Q T u>. 
In the base state, the rod is straight with r = zz, the 
tension A vanishes, and the rod undergoes rigid-body 
rotation about the z-axis with speed lo = ujq, leading to 
a constant gradient in twist. At the free end, s = L, the 
torque vanishes: Cil = 0. Thus, 



(±8tT± = -Av"f + Czx (Qr'i)' . 



(72) 



n = Crw (s - L)/C. 



(70) 



To estimate when the distributed viscous drag causes 
the rod to twist so much that the straight state be- 
comes unstable, we may consider the related problem 
of a rod twisted by external moments acting at either 
end, in the absence of viscous forces. In this simpler 
problem, the straight state of the rod becomes unsta- 
ble when the twist moment is roughly equal to the 
characteristic bending moment of a rod of length L: 
Cf2 crit ~ A/L (lLandau and Lifshitzl. 1 19861). S ince C « A 
for ordinary materials ( Landau and Lifshitzl . fl986l ). the 



where rj_ = (x, y, 0) and the twist is given by Eq. ([701) . To 
leading order, the deflection of the rod is perpendicular 
to the z-axis. We therefore expect that the tension A 
is second order in x and y, which may be confirmed by 
expanding Eq. (|66p in powers of x and y. 

Perturbations to the twist are governed by the diffu- 
sion equation with diffusion constant C/£ r . Thus, the 
perturbations in shape and twist may be studied inde- 
pendently. Since perturbations to twist are not destabi- 
lizing, we focus on Eq. (|72p . which governs the shape of 
the filament. As in the static twist instability, this equa- 
tion is most conv eniently analyzed using com plex nota- 
tion £ = x + iy ( Landau and Lifshitzl . [1986). In these 
variables, Eq. (f?2"|) becomes 



CLfit£ = -^"" + i&wo [(s-L)?]'. 



(73) 



For boundary conditions, we assume the end at z = is 
held fixed, £(0,i) = 0, and clamped, £'(0,i) = 0. At the 
far end the conditions are vanishing force, £/"(L,t) = 0, 
and vanishing bending moment £"(L, t) = 0. The so- 
lutions to Eq. (|73p are proportional to exp(icrf). When 
the imaginary part of a is negative, the base state is 
unstable; the real part of a is the rate of rigid body ro- 
tation of the filament about the z axis. Equation (|T3"1) 
has nontrivial solutions only for special values of loq; to 
determine the value for which the straight state first be- 
comes unstable, we assume £(s,t) = expi^i, where x 
is real, and use the boundary conditions to determine 
the critical uip. This pro cedure is best done numerically. 
IWolgemuth et all (|2000f) found 



Wcrit « 8.9 



A 



(74) 



in agreement with our rough estimate Eq. (|71[) . The mo- 
tion of the filament can be described in terms of two dif- 
ferent rotations. One is the rapid local rotation of each 
element of the rod about the local tangent vector; this 
rotation has rate loq. But the center line of the filament 
undergoes rigid-body motion at a different, much smaller 
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rate Xt which at u> — w C rit takes the value 2 

.4 



X «22.9 



(75) 



The rigid-body rotation rate x is °f the same order of 
magnitude as the relaxation time for bending. The ratio 
of cj cr it to x is °f order L 2 /a 2 since the ratio of the corre- 
sponding resistance coefficients Cl/Ci- is of order L 2 /a 2 . 

Numerical methods applied to the full equations of mo- 
tion are required to find the whirling state for uiq > w cr ;t- 
IWolgemuth et all ( 2000f) attempted a weakly nonlinear 
ana lysis by applying a num erical pseudospectral analy- 
sis ( Goldstein et al. , [l996a) to the equations of motion 
truncated at second order in displacement, and found 
a small-amplitude steady-state whirling state for uj just 
above w cr it- However, subs equent work using the im- 
mersed boundary method ( Lim and Peskinl |2004 ). as 
well as simulations that include the full coupling of ther- 
mal, elastic, and hydrodynamic forces ( Wada and Netzl . 
I2006D . have shown that there is no small-amplitude 
whirling state, and that the weakly nonl i near a nalysis 
yields incorrect results. IWada and Netd ( 20061 ) found 
the s ame value fo r the c ritical frequency that we quote 
here: IWokemuthl (I2009D has also found this value as 
well as the large amplitude whirlin g state using the 
pseudospectral method s empl oyed in ( Wolgemuth et all , 
I2000D . [L im and Peski n| (|2004l ) found a much lower value 
for the critical frequency. The reasons for the discrepancy 
remain unclear. 



IV. MEMBRANES 

In this section we follow the same format as our dis 
cussion of cu r ve dy namics, but we consider fluid mem- 
branes (jBoall , 120021) . 



We begin with the geometry of a 
fixed surface in space. This material is standard, but we 
include it for completeness. Then we turn to the cal- 
culation of membrane forces per unit area as a function 
of membrane shape. Since traditional approaches to this 
calculation are somewhat tedious, we revie w a recent ide a 
of Guven that embodies economy of effort (|Guvenl . l2004r ). 
We then revisit the basic geometric concepts for a moving 
surface, taking care to distinguish between time depen- 
dence arising from flow or surface deformation and time 
dependence arising from time-dependent coordinates. Fi- 
nally we conclude with the equations of motion for an 
incompressible fluid membrane and apply them to the 
instability of a cylindrical membrane tube under pres- 
sure. 



A. Geometry of a stationary surface 

The intrinsic geometry of a space curve is Euclidean: a 
one-dimensional being confined to a space curve cannot 
tell if the curve bends and twists in the ambient three- 
dimensional space. This fact is why we can always choose 
arclength coordinates for a curve. In contrast, curved 
surfaces are non-Euclidean. Anyone who has attempted 
to flatten an orange peel onto a tabletop knows that it is 
impossible to construct a Cartesian coordinate system on 
a sphere. Thus we are forced to use curvilinear coordi- 
nates when a surface has curvature, and the description 
of the geometry of surfaces is necessarily more compli- 
cated than for curves. A key theme is that geometrical 
(and physical) concepts do not depend on the choice of 
coordinate system. Thus, we will gradually be led to the 
basic ideas of tensor calculus, which will allow us to con- 
struct geometrical quantities such as curvature in terms 
of any curvilinear coordinate system. 

An additional complication arises for fluid surfaces 
that was not present for solid rods: even for a stationary 
shape, there may be a complex flow pattern on the sur- 
face. We will return to this complication later; first we 
will consider stationary surfaces and review the geome- 
try necessary to calculate the force per unit area due to 
internal stresses as a function of membrane shape. 

Denote the coordinates of a stationary surface as 
i? 1 ;^ 2 }- Greek letters a, f3, ... refer to the two- 
dimensional coordinates, such as £ Q , and Latin letters 
i, j, ... continue to refer to Cartesian three-dimensional 
coordinates. Thus the position of a point in space on 
the surface is given by r^ 1 ,^ 2 ). The distance between 
two points on the surface is given by the Pythagorean 
theorem 



ds 2 = dr 2 = -^- ^dCd^. 



(76) 



The right-hand side of (jT6"|) is known as the first funda- 
mental form. The quantity 



= d a r ■ d$r 



(77) 



2 IWolgemuthl (I2009T) has provided a more accurate value than the 
one stated in j Wolgemuth et al.l . l200oh . 



is known as the metric tensor, because it converts dis- 
tances {d^,d^ 2 } measured in coordinate space to dis- 
tance measured in real, three-dimensional space. Since 
distances can be measured by a being on a two- 
dimensional surface, the metric tensor is an intrinsic 
quantity. 

A surface has two independent tangent vectors at any 
point, which may be taken to be t Q = d a r along the 
direction of increasing £ Q , for a = 1,2. In general, the 
basis {ti,t2} is not an orthonormal basis. The tangent 
vectors t a are not geometric objects, since they depend 
on the choice of coordinates, but a tangent vector field 
W = W a t a is a geometric object independent of the 
choice of coordinates. The unit normal to the surface is 
given by the cross product of the two coordinate tangent 



14 



vectors, suitably normalized: 

ti x t 2 



n 



ti X to 



(78) 



Note that there are two unit normal vectors at every 
point, n and — n. Unless there is an asymmetry such 
as different fluids on either side of the membrane, there 
is no reason to prefer n over — n. The element of area 
on a surface is given by the parallelogram formed by the 
infinitesimal vectors txd^ 1 and t 2 d£ 2 , which has area 



dS = n- (t^e 1 x t 2 d£ 2 ) = Vsd^ 2 



(79) 



where g — gwg-n.— 512921 is the determinant of the metric 
tensor. 

The curvature of a space curve is the tangential com- 
ponent of the rate of change of the normal to the curve 
[Eq. (|22p]. Since a surface has two tangential directions, 
curvature of a surface is described by a bilinear form 
known as the second fundamental form: 



dn ■ dr = -K af3 dCd^ 



(80) 



where 



K, 



a/3 



-d a h ■ t 



fi- 



ll) 



Note that since n • tp = 0, we have 

9 Q n • tp = — n • d a tp = — n • dpt a = dph ■ t a , (82) 

and therefore conclude that the second fundamental form 
is symmetric: K a p — Kp a . 

To make a stronger connection between our discussion 
of curves and the second fundamental form, consider a 
curve on a surface, parameterized by arclength, x(s) — 
r(^ 1 (s), £ 2 (s)), and define the normal curvature K n by 
K n = n • kN = fi • dT/ds, or 



d£ Q d fd£_ 
~ds~dt[ 



cfi 



ds 



d^d^ 

= ^ IT dT" 



(83) 
(84) 



The normal curvature is the ratio of t he second and the 
first fundamental forms (jStruikl . [l988l ). 



K a pd^de 
<7^d£Me ' 



(85) 



Equation (|85[) implies that at a given point, all curves 
with the same direction d^/d^ 1 and parameterized by 
arclength have the same normal curvature. 

It is useful to characterize the curvature of a surface 
by the maximum and minimum values of the normal cur- 
vature at a point. If v = v a t a is a unit vector, then the 
extreme values of K n are given by the critical points of 



/ = v a K afjV P -\{v a g a ^ -1), 



(86) 



where A is a Lagrange multiplier. 
df/dv a = yields 



K af3 v p 







Demanding that 
(87) 



Note that multiplying Eq. (|87j) by v a and summing over 
a, together with the definition of normal curvature (|85[) . 
implies that A = K n . Since v a ^ 0, the condition for 
Eq. (|87|) to have a solution is det(K a p — n n g a p) = 0. If 
we use g a @ to represent the inverse of the metric tensor, 



9 ai 9~/fi 



(88) 



then what we have shown is that the maximum and min- 
imum values of the normal curvature are given by the 
eigenvalues of the matrix 



K a g = g^K lf3 . 



(89) 



These extreme values are known as the principal curva- 
tures of the surface. Not that since K a p and g a p are 
symmetric, the principal directions v a corresponding to 
the two principal curvatures are orthogonal. Note that 
Eq. (f8Tj) may be written in terms of K a ^ to yield the 
Weingarten equations: 



(90) 



As an aside, we emphasize the importance of whether 
or not an index is raised or lowered. The convention is 
to use the metric tensor to lower indices, and the inverse 
of the metric tensor to raise indices. Quantities with 
upper and lower indices transform oppositely from one 
another under coordinate changes; for example, W a and 
t Q transform under under coordinate changes in such as 
way that the vector field W = W a t a remains invariant. 
These transformation rules are reviewed in Appendix [A"l 

Returning now to the discussion of curvature, we note 
that it is common to describe the curvature of a surface in 
terms of the mean curvature, H = K a a /2, and the Gaus- 
sian curvature, K — detK a ^. The mean curvature H is 
the average of the sum of the two principal curvatures of 
a surface. It is a measure of how a surface bends in three- 
dimensional space; H is analogous to the curvature k of 
a space curve. A plane has vanishing mean curvature, 
a cylinder has mean curvature equal to half the inverse 
of its radius, and a sphere has mean curvature equal to 
the inverse of its radius. The Gauss curvature K is the 
product of the two principal curvatures of a surface. It 
is an intrinsic quantity; K can be determined by two- 
dimensional beings making local measurements of length 
on the surface. For example, by choosing appropriate co- 
ordinates, one can show that the relation between the cir- 
cumference C a nd radius R ofa s mall circle on a surface 
involves K only (|do Carmol[l976l ): C « 27ri?(l- KB 2 /6) 
(Fig. [1]). Below we use direct calculation to see that K 
is an intrinsic quantity. Note that a plane and a cylinder 
both have vanishing Gauss curvature, which is consis- 
tent with the fact that local measurements confined to 
the surface cannot distinguish a plane from a cylinder. 
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FIG. 8 (Color online.) The height near the origin is given 
by the second fundamental form: h ~ K a pA^ a A^, where 
A^ 1 = Ax and A£ 2 = Ay. 



In our study of curves we found it useful to express the 
derivative of directors in terms of the directors. Likewise, 
it is useful to examine the derivative of the coordinate 
tangent vector fields t a in terms of t a and n. Since the 
coordinate tangent vectors on a surface need not have 
unit length, the derivative of a tangent vector has both 
tangential and normal components: 



d a t 



K, 



a/3 1 



r 7 t 



(91) 



where the normal components have been identified using 
the Weingarten equations (|9H)) and n • t Q = 0, and the 
as yet unknown quantities r 7 ^ are called the Christoffel 
symbols. Note that r 7 ^ = TZ a since d a dpr = dpdar and 
K a fj — Kp a . The formula K a p = n • d a dpv also suggests 
another interpretation of the second fundamental form. 
Choose a point p on the surface and introduce Cartesian 
coordinates {x, y, z} in three-dimensional space such that 
z is parallel to n. Near p, we may choose coordinates 
{C 1 ;^ 2 } ~ { x iV}- The height of the surface over the 
tangent plane is given by the second fundamental form 
(Fig.©. 

To see how the Christoffel symbols are related to the 
shape, observe that 



da9pj = ^(t^ • t 7 ) 



r a /35<57 



(92) 



By cyclically permuting a, (3, and 7 in Eq. (|92[) and then 
adding and subtracting appropriately, we can solve for 
the Christoffel symbols to find 



1 aB 



g lS r 5al} 



(d a g~fp + dpg^a - d^gafj) . (93) 



Since the Christoffel symbols depend only on the metric 
and its derivatives, they are intrinsically quantities. 

The equations (|9 1 [) are sometimes called the Gauss 
equati ons or the p artial differential equations of sur- 
faces (Struik, 1988). These equations are indispensable 
for formulating the equations of motion for fluids on a 



deforming surface. In particular, we will need to take 
derivatives of the velocity field, which has tangential and 
normal components. However, we will limit this section 
to a simplified discussion in which we consider the deriva- 
tive of a purely tangential vector field, W = W a t a : 



8 a W = (d a WP) tp 



^r 7 ,t 7 



W K a 



/?n. 



(94) 



Once the direction t a of the derivative is chosen, the tan- 
gential and normal components of <9 Q W are uniquely de- 
fined, independent of choice of coordinates. These quan- 
tities contain all the information about how W changes 
along the direction t^. In contrast, the quantities daW 13 
do not contain complete information about the change in 
W. For example, we can imagine a vector field on small 
region of a curved surface in which the components W 13 
are constants. Then d a W^ = 0, but W is not constant, 
since the surface is curved. Therefore, we are lead to de- 
fine the covariant derivative of the vector field W@ as the 
tangential components of 9 Q W: 



(y a wf = d a we + r^w. 



(95) 



The covariant derivative of a tangent vector field is an 
intrinsic quantity since the Christoffel symbols are in- 
trinsic. 

The covariant derivative is readily generalized to other 
quantities defined on a surface. For scalar functions, 
V Q / = d a f. Since each Cartesian component of a vector 
field may be considered a scalar function, we may also 
write V Q W = 9 Q W. The Leibnitz rule is used to define 
the covariant derivative of higher-order tensors. For ex- 
ample, we may write the dot product of two vector fields 
U and W as U • W = U a W a , where U a = g^U 13 . The 
Leibnitz rule then implies 



V a (U p W ) = (VaU)pW^ + U p {V a Wf. 

Since UpW 13 is a scalar function, V a ([/^W /3 ) 
daiUpW 13 ) and solving © for (V Q W> yields 

Since the Cartesian coordinates of a vector are scalar 
functions, we may write 



(V a W) = d a W 



(96) 



(97) 



r I / 3 t 7 = K apri- 



(98) 



At first sight this relation may seem confusing since we 
defined the covariant derivative of W 13 by removing the 
normal part from d a W, and here we have normal compo- 
nents in the covariant derivative of tg. The key point is 
that the two-dimensional indices determine the nature of 
the covariant derivative. The rule (|9"8]) ensures that the 
Leibnitz rule holds for the components and basis vectors 
of a vector: V Q W = d a W = (V a VF)% + W p V a tp. 

The general rule is that the covariant derivative of a 
tensor has a term with a Christoffel symbol added for 
every upper index, and a term with a Christoffel symbol 
subtracted for every lower index. For example, 



r k 5 



(99) 
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Since the metric tensor g a p is derived from the Euclidean 
metric of the ambient three-dimensional space, the co- 
variant derivative of the metric tensor vanishes: 



V, 



= d r 



■9M 



O/w = o. 



(100) 



This fact may be deduced by writing Eq. ()96|) in terms 
of U a , W 13 , and g a p. The metric tensor and its inverse 
is free to pass through the covariant derivative. 

To formulate the generalization of the Stokes equations 
(|57H58[) for membrane flow on a curved surface we will 
need to consider multiple covariant derivatives acting on 
the velocity field. Special care must be taken, since co- 
variant derivatives do not commute when acting on a vec- 
tor field. The Gauss-Codazzi equations, which we now 
derive, determine the relation between (VqV^)!^ 7 and 
(VpV a )W J . Since the tangent vector field W is three 
scalar functions on the surface, covariant derivatives com- 
mute when acting on W: 



(V a V /3 - V V Q )W = 0. 



(101) 



Resolving this formula into tangential and normal com- 
ponents yields the Gauss-Codazzi equations: 

(V Q V^-V^V a )^ = {K\K p S ~K\K aS )MiQ2) 
V a K Pl = VpK aj . (103) 

The right-hand side of Eq. (jl02|) can be simplified by 
noting that K^ a Kp —K J ^Kj is asymmetric in (a, (3) and 
(7,(5), and therefore proportional to the determinant of 
K°g. The only independent component is a = 1, /3 = 2, 
and we may write 

(ViV 2 - V 2 Vi)iy 1 = KW 2 (104) 
(V1V2 - V 2 Vi)VK 2 = -KWi. (105) 

The Gaussian curvature, originally defined in terms of 
the extrinsic quantity K a p — n • dadpr, is intrinsic. The 
fact that K may be determined by two-dimensional be- 
ing living on a surface i s known as Gauss's Theorema 
Egregium (jMorganl . 119981 ) . 

We close this subsection with a discussion of the co- 
variant form of Green's theorem for a two-dimensional 
surface. This relation will be required when we derive the 
stress tensor for a membrane. First note that Eq. f|93[) im- 
plies that the covariant divergence of a two-dimensional 
vector field has a simple form: 



\9 



(106) 



This formula is useful when integrating by parts, since 
the area element dS = *Jgd 2 t; depends on £ a : 



V a d a fdS = 



J(V a V a )fdS + J V a (V a f)dS. 

(107) 



In appendix [Bl we show that the second term of Eq. (|107[) 
can be written as an integral over the boundary of the 
surface: 



W a B a dS 



B a Pa ds, 



(108) 



i)D 



where p = dr/ds x n is the unit tangent vector perpen- 
dicular to the boundary 3D. 



B. Energies and forces 

In this section we find the force exerted on a patch 
of a membrane by the rest of the membrane, or, equiv- 
alently, the force the patch exerts on the surround- 
ing medium. We follow the geometric approach o f 
IGuvenl (12004): see also (ICapqvilla and Guvenl l2002f ) 



and ( Capovilla and Guvenl . 2004h . For concreteness. we 
focus on fluid membranes; however, the approach may 
also be applied to memb ranes with internal degrees of 
freedom such as tilt order ( Miiller el aZ.I . l2005h . The elas- 
tic energy for a fluid membrane takes the form 



E 



I [3 < 2ff > s 



RK 



AS, 



(109) 



wher e k and k are mo duli with units of energy ( Canhaml . 
119701 : iHelfrichl . Il973l ). For a closed vesicle, the Gauss- 
Bonnet theorem implies that the integral of the Gaus- 
sian curTOture_K_over the surface is a topological invari- 
ant ( Kamienl . [2002) . and therefore does not change under 
a variation of the shape. Henceforth we drop this term. 

To find the force f dS on a patch for an arbitrary shape, 
we invoke the principle of virtual work as we did with 
filaments. We suppose the shape is held in equilibrium 
by an external force per area — f(^ 1 ,^ 2 ) and examine the 
condition SE to t = 0, where 



SE, 



SE 



f • SrdS. 



(110) 



Finding the variation of E for surfaces is harder then 
it is for curves. Since E = J £(g a/ 3, K a/ 3)dS , we must 
find 5g a p and 5K a p in terms of Sr. This direct approach 
is str aig htforwa r d but tedious (|Zhong-can and Helfrichl . 
11989b . IGuvenl {2004) suggested introducing Lagrange 
multiplier functions which allow the variations r, g a p : 
and K a f3 to be taken independently: 



SW = SE + Jf- SrdS + S{ J F a ■ (d a r - t a )dS 
+ J [Ait ct -n + A n (n 2 -l)]d,5' 

- (T al3 /2)(g aP -t a -t )]dS) . 



Ill 



The energy density £ is now regarded not as function 
of r but instead as a function of g a p and K a p. The 
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condition of equilibrium is that SW = for the variations 
5r, St a , r)n, 5K a p, and 5g a p. Since the variables appear 
quadratically in the constraints, the variations are easy 
to compute. For example, the variation with respect to 
r implies that 



V a F a = f . 



(112) 



Thus the Lagrange multiplier function F a is the two di- 
mensional stress tensor. To see why the stress is of mixed 
character, it is useful to apply Green's theorem (|108p : 



V Q F Q dS = * F a Pa ds. 

D JdD 



(113) 



The force per unit length acting through a contour in 
the surface is a three-dimensional vector, and the normal 
p to the contour lies in the local tangent plane and is 
therefore a two-dimensional vector. 

It is convenient to define the functional derivative 
SW/Sh a , where h Q is a generic mixed tensor, as 



SW 

SW= I — ■ Sh a dS. 
Sh a 



(114) 



The condition SW/6t a — yields the stress in terms of 
Lagrange multipliers, 

F Q = (T a!3 - A aj K^) tp + Xfn. (115) 

Likewise, the condition SW/Sri = along with Eq. 
gives A" and A n in terms of A Q ^: 



XI = VpA a P 
2A n = K aP A a ?. 

Finally, 8W/SK a p = implies 

A a/3 = _ £a 



(116) 
(117) 



(118) 



where £ a0 = d£/dK a0 . 

We need a few identities to help us find the variation of 
E with respect to the metric tensor. For example, to find 
the variation of the determinant of the metric tensor, we 
use the identity 



exp log det g a p = exp tr log g a p, (119) 



which yields 



5g = gg aP 6g a(3 . (120) 



Thus SdS = g af3 Sg a pdS/2. Note also that Eq. (g5} im- 
plies 

Sga p = -g^Sg^sg 60 . (121) 
With these formulas we see that SW/ Sg a p — implies 

0£ 



T a/3 _ g a0 £ + 2 



dg a 



122) 



The two-dimensional stress tensor may also be written as 



d 



(123) 



Note that the contributions from varying dS in the con- 
straint terms vanish since the constraints are invoked af- 
ter taking the variation. To summarize, we have found 
a general formula for the stress tensor of a membrane in 
terms oiE a ^ = d£/dK a p, the cu rvature tenso r K a p, and 
the two-dimensional stress T a @ (jGuvenl . l2004h : 



(T 



a/3 



(Vp£ a ?)h. 



(124) 



The variational approach with auxiliary variables can 
also be used to find the moments per length act- 
ing ac ross any infinitesimal line segment in the mem- 
brane (|Capovilla and Guvenl . l2002t iMuller et all , l2007fh 
Two examples illustrate the formalism. 



1. Soap film 

The energy of a soap film is given by the interfacial 
tension 7 times the area. Thus, the energy density £ — 
7, and £ al3 — 0. The stress depends only on the two- 
dimensional tensor T a @ = jg 01 ^ , and 



film 



19 ap t p . 



(125) 



Fg lm is in the tangent direction as expected. The force 
per unit area acting on the film is ffi ml = V Q Fg lm or 



Ifiln 



2 1 Hh+[d al )g aP t 



0- 



(126) 



The tangential component of Eq. (|126|) corresponds to a 
Marangoni stress and arises when surfactant concentra- 
tion varies over the film. For soap film on a frame in 
equilibrium, 7 is uniform and f = is the minimal sur- 
face equation H = 0. If the film forms a closed surface 
with fixed enclosed volume, we introduce another La- 
grange multiplier as in Eq. (|27[) and find that the mean 
curvature must be constant. 



2. Fluid membrane 



Consider the fluid membrane energy (|109[) with k = 0. 
The energy density is £ = (n/2)(g a l 3 K a p) 2 . To find 
we must differentiate £ with respect to g a p while holding 
K a p fixed. Note that Eq. (|121[) implies 



dg 



(127) 







where we have used the symmetry g a p = gp a . A short 
calculation reveals that the bending stress in a fluid mem- 
brane is 

Fg ond = 2k [(g^H 2 HK^) t p - (W a H) n] . (128) 
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The force per unit area acting on the membrane due to 
internal bending stresses is given by the divergence of the 



stress tensor F^cnd' 



f bcnd = -2k (V 2 H + 2H 3 - 2HK) n. 



(129) 



We used Eq. (|103[) to show that the tangential compo- 
nents of fbond vanish for constant k. There is no tan- 
gential component of the force per unit area because the 
bending energy is invariant under changes in coordinates, 
and a small change in coordinates corresponds to a defor- 
mation of the surface along tangent directions at every 
point, 



(130) 



Equation (|129jl for fbcnd may be derived directly 
from f = —SE/dr with con siderably more ef- 
fort (|Zhong-can and HelfricfJ . Il989h . 

Although we argued that membranes in the tense 
regime are approximately incompressible, it is instruc- 
tive to examine the stress and force per unit area aris- 
ing from the stretching energy, E s — J £ s dS, where 
£ s = (fc/2)(0/ '<f>o — l) 2 . First we will work by anal- 
ogy with a two-dimensional liquid and identify the two- 
dimensional pressur e P2d by considering a situ ation with 
uniform density <j> (|Cai and Lubenskvl [l995). In that 
case, the two-dimensional pressure is p2d = —dE s /dA, 
where A is the area and the derivative is taken at fixed 
number <pA. Carrying out the differentiation leads to 
P2d = — (£s ~ 0£s)j where £ s is regarded as a function of 

To apply the principle of virtual work to find the stress 
and elastic force, we must keep the identity of the par- 
ticles fixed under the variation. This constraint arises 
because the principle of virtual work is derived by as- 
suming mechanical equilibrium of a fixed set of particles. 
Therefore, the variation satisfies 8(yJ~g<j)) = 0, which im- 
plies 5(f> = -{l/2)g af3 8g af3 (t> and 

6E S = - y 1 y/gg ap 5g a0PM dS. (131) 

From Eq. (|122[> . the corresponding two-dimensional 
stress tensor is T"^ = —p2dg al3 and we see that the two- 
dimensional pressure acts like a negative interfacial ten- 
sion [compare with Eq. JTSSJl ]. Thus, F S Q = -_p 2 d5 Q/3 t /3 
and f s = -2p 2 dHh - (d a p 2 d)g al3 t- 

We can make a further simplification by observing that 
the large modulus k implies that the density will be close 
to its preferred value: <j> — 4>o + <f>i, where 4>\ <C <po- 
Working to first order in </>i leads to p2d = &<^i/<fo- The 
stress arising from stretching is therefore 



F? = -(k<h/<Mg a % = 
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t/3, 



(132) 



where 7 = —p2d- In general, the stress F" cannot 
be disregarded in a membrane, since the smallncss of 
the stretch <f>i is offset by the largeness of the modulus 
k. Note that in mechanical equilibrium, the tangential 




FIG. 9 Equilibriu m shapes of a vesicle su bject deformed with 
an optical tweezer ijFygenson et al.Lll997l h The magnitude of 
the force increases from left to right. Figure courtesy of D. 
Fygenson. 



stresses acting on the membrane are constant, and there- 
fore p2d and 4>\ are uniform, independent of the shape of 
the membrane. Often, membrane shapes are studied un- 
der the assumption of incompressibility, k — > 00, in which 
case 7 = — P2d is treated as a constant Lagrange multi- 
plier, and V Q F" is the associated force per area area 
enforcing the constraint. For dynamics, the Lagrange 
multiplier may have spatial dependence if there are tan- 
gential stresses. 

The formulas for the stress tensors Fg d (|128p and 
F" (|132|) give insight into the mechanics of membranes. 
Consider a tu bule pulled out of a ves i cle with a point 
force , (Fig. |S) (E vans and Yeunsd . Il994t iFygenson et all 
I1997D . It is common to regard the spherical part of the 
vesicle in the right- most frame of Fig. [5] as a reservoir of 
lipid, and to suppose there is an interfacial tension 7 that 
remains constant during the deformation. Equivalently, 
we may regard the membrane as under tension, with 7 — 
~P2d > and constant. Since the spherical part of the 
vesicle is much bigger than the tubule, we suppose that 
P2d is independent of tubule length once the tubule is 
formed. Thus, £ = 7 + {n/2){2H) 2 . The cylindrical 
part of the vesicle is easy to analyze, since K = and 
H = — l/(2R), where R is the tubule radius. We may 
disregard the pressure jump across the membrane since 
the large radius of the spherical part of the vesicle ensures 
that the pressure jump is small. The total force on an 
infinitesimal patch of the cylinder therefore vanishes, f — 

ffilm + fbend = 0, Or 



j/R+k/{2R 3 ) = 0. 



(133) 



In equilibrium, R = ^k/(2j). 

Closer inspection of the equilibrium configuration may 
seem to reveal a paradox dFournier , l2007t iPowers et all 
120021 : IWaueh and Hochmuthl . 1198711 . With the usual 
cylindrical coordinates {f 1 ,^ 2 } = {<£, z}, we expect that 
the interfacial tension 7 will lead to a force per unit 
length F v = -f(p along the edge AB of Fig. [101 with a 
force per unit length of equal magnitude along the other 
edge, CD. But since there is no pressure jump across 
the surface, there are no other forces to balance these 
forces. To resolve this apparent paradox, examine the 
stress F a = Fg lm + F bond . In cylindrical coordinates, the 
first fundamental form is ds 2 = R 2 d(p 2 +dz 2 , and the sec- 



ond fundamental form is K a pdl; a dl; 



-Rdtp 2 . Along 
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FIG. 10 (Color online.) A section of a tubule, showing the 
forces per unit length F z and F 1 " 3 . 



a circular contour of constant z, F z — 2jz. Thus, the 
force required to hold the tether in equilibrium is Attjz; 
twice what it would be for a cylinder with no bending 
stiffness. However, along a line of longitude, constant 
ip, F v = [7 - k/(2R 2 )]<p = 0; the contributions to the 
stress from the bending energy cancel the contributions 
from the interfacial energy, as they must since there is no 
pressure difference. 



C. Moving surfaces: kinematics 

1. Coordinates on deforming surfaces with flow 

The description of a moving surface which may have 
internal flows requires careful thought. Curves are much 
simpler. We saw that we can always use arclength for the 
spatial coordinate, and in the case of an inextensible fil- 
ament, the arclength coordinate of a material point does 
not change. There is no analog of arclength parameter- 
ization for surfaces, and instead we must choose coordi- 
nates that are most appropriate for the situation. There 
are several common choices, including convective or ma- 
terial coordinates, surface-fixed coordinates, and general 
coordinates. Material and surface-fixed co ordinates ar e 
useful for formulating equations of motion ( Arid . fl989t) . 
while general coordinates are useful when one must solve 
for the evolution of a shape. In this section we briefly 
describe the most important aspects of each class of co- 
ordinates. 

To define material or convective coordinates, 
parametrize the surface at time t = with coordi- 
nates u r = {u 1 ,!! 11 }. We will use Roman numerals or 
upper-case Greek letters for material coordinates. Since 
we disregard thermal effects such as Brownian motion of 
the molecules making up the membrane, the motion of 
the material points defines a smooth flow velocity. The 
coordinates u r are labels that stay with the material 
particles as time evolves, with the trajectory of the 



material point {u^u 11 } given by r(w / ,u // ,i). The 
velocity of a material point is given by 



As mentioned above, convected coordinates are useful for 
formulating the equations of motion. For example, we 
will see in flIVD that the Reynolds transport theorem 
for conservation of parti cles i s conv eniently derived using 
convective coordinates (lAris]. H989h . However, since the 
motion of a surface is typically something to be solved 
for, r(u ! ,u n ,t) is usually unknown at the outset of any 
problem. 

An alternate approach is to attempt to define a coor- 
dinate system that is fixed to the surface, independent of 
the flows of the molecules making up the surface. Such a 
choice is easy enough for a surface that does not change 
shape, but is impossible for a deforming surface. The best 
we can do is to choose coordinates £ 2 } such that the 
velocity of a point with fixed coordinate — not necessarily 
a material point — is purely normal to the surface: 



dr 

di 



= 0. 



Thus, 



V = 



dr 

dt 



= t a V a + Kn, 



(135) 



(136) 



where V a = d£ a /di, £, a (t) is the trajectory of a par- 
ticle with constant u r , and V n h — dr/dt. Note that 
dr/dt is purely normal. Surf ace- fixed coordinates are 
conceptually useful when formulating equations of mo- 
tion of moving surfaces. For example, we will see below 
how the condition Q135P leads to alternate expression the 
contin uity equation that commonly appears in the liter- 
ature (lAridJl989l:lBuzzal.l2002l : lDorries and Foitm1 . il 9961: 
iFuiitanil . Il997t lYouheil ." Il994f ). But like material coordi- 
nates, surface-fixed coordinates also require knowledge of 
the shape evolution which is usually unknown until the 
complete problem is solved. 

Therefore, it is common to use general coordi- 
nates when one is solving for the shape of a mem- 
brane as a function of time. By general coordi- 
nates we mean coordinates in th e Eulerian point of 
view ( Chaikin and Lubenskvl . Il995f) , which are not con- 
strained by conditions such as (|135f> or the Lagrangian 
label of convected coordinates. For example, we will use 
cylindrical coordinates r(6, z,t) — h(z, t)r + zz in our ex- 
ample of axisymmetric surfaces. (Here f is a unit vector 
in the radial direction in cylindrical coordinates.) The 
velocity in general coordinates is 



dt 



dr du a dr 
du"~dT + W 



(137) 



where dr/dt has both normal and tangential components. 



2. Strain rate 



v= - 

dt 



= v T t T 



Kn. 



(134) 



To formulate the hydrodynamics of fluid membranes, 
we need to measure the strain rate, defined as rate that 
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the distance between two nearby points on the surface 
changes with time. If dr is the vector connecting the two 
nearby points at a given time, then after a time interval 
dt, the two points are connected by dr' = dr + (dV)di. 
The strain rate tensor S a p is then defined by 



ds' 2 = ds 2 + 2dtd^ a d^S a ^ 



(138) 



where ds' 2 = dr' • dr' and ds 2 = dr • dr. Writing ds' 2 in 
terms of dr and V leads to 

S a p = \ {d a V ■ tp + d fj V • t ) . (139) 

Note that S a p is unchanged by rotation with constant, 
uniform r ate u), V = us x r, as requ ired by the definition 
of strain (jDorries and Foltirl 1 199(a ). Likewise, the strain 
is invariant under translations r i— ► r + a and Galilean 
boosts V i— > V + Vn , where a an d Vo are constant and 
uniform ( Dorries and Foltirl Il996f ) . Applying the Gauss- 
Weingarten equations (|91l90p to Eq. (|139[) yields an al- 
ternate form, 



S a p = ^ (V«^s + VpV a - 2K a pV n ) . 



(140) 



The first two terms the right-hand side of Eq. (I140j) 
are the covariant generalization of the strain rate of a 
two-dimensional material. The third term represents the 
strain that takes place for example when a sphere inflates 
uniformly — even though there are no velocity gradients, 
there is strain whenever there is curvature and a normal 
velocity V n . Note that gradients in V n do not lead to 
strain at first order in dt. 



3. Equation of continuity 

It is useful to formulate the covariant version of the 
conservation of particles, or, more generally, formulate 
the covariant v ersion of the Reynolds transport theo- 
rem ( Arid . [l989l) . To this end, consider a patch of mate- 
rial points S(t) on a fluid membrane, and a scalar quan- 
tity <f>, which for concreteness can be taken as the number 
per unit area. We wish to calculate the rate of change of 
the total number of particles in the patch S(t), assuming 
that the particles never leave the membrane to enter the 
solvent. To eliminate the time-dependence of the patch 
region, use convected coordinates: 



iL [ <£dS=£/ ^{vF^VG&iSdu 11 , (141) 
M Js(i) dt J Sg 

where Sq = S(t = 0), and G is the determinant of the 
metric in the material coordinates {u^u 11 }. The total 
time derivative of <f> is given by the chain rule, 



dt 



di 



V a d a 



(142) 



where V a is the tangential component of the veloc- 
ity V = V a t a + V n h. To evaluate the time deriva- 
tive of VG, note that Eq. (fT2"Tj) implies dVG/dt = 



VG(G rA dGrA/dt)/2, and dGr A /dt = <9 r V ■ t A + t r 
9a V. Reverting to general coordinates, we have 



G 



TA 



dG rA 
dt 



n a/3 



(d a V ■ tp + dp V • t a ) 



(143) 



We may now use the Gauss- Weingarten equations (|91I90|) 
to deduce 



d 

dt 



<t>dS 



S(t) 



S(t) 



dt 



+ <jN a V a - 2(j>HV n 



-±+(j)V a V a 
dt 



2<pHV a = 0. 



dS, 
(144) 

(145) 



The continuity equation (|145[) is sometimes expressed 
in terms of derivatives o f the metri c tensor using 
surface-fixed coordinates 



Buzza 



iDorries and Foltinl . Il996t iFuiitanil . Il997t lYouheil 11994 
in which the condition (|135|) can be used to show 



j_d_ 

Vgdt 



(146) 



However, since Eq. (1 146[) only holds for surface-fixed co- 
ordinates, and not for example in the cylindrical coordi- 
nates we use to solve our illustrative problem in ^IV.El 
we will use Eq. (|145[) in this review. 

For an incompressible membrane, <j) is constant, and 
the continuity equation simplifies to 



V n V a = 2HV„ 



(147) 



In terms of the strain rate, incompressibility implies 
S« = 0. 



D. Dynamical equations for an incompressible membrane 

There are several complementary approaches to the 
dynamics of membranes from the points of view of 
fluid mechanics and interface s d Arroyo and DeSimone , 

19601: IWaxman 



20091: iKralchevskv et all 11994 IScrivenl 



1984 ) . statistical mech anics (|Cai and Lubenskvl . 1994 , 



1995; Mi ao et all [2002] ) . and the director approach anal- 
ogous to our treatment of rods ( Hu et al\ , [2007) . In- 



stead of presenting the most general models, we dis- 
cuss the simple case of an incompressible fluid mem- 
brane. We disregard the bilayer nature of the membrane 
and t herefore do not consider effects such as bilayer fric - 
tion (jSeifert and Langerl Il993b lYeung and Evan s, 1995). 

The form of the two-dimensional viscous stress ten- 
sor is derived the same way as the stre ss for three- 
dimensional fluid mechanics in flat space ( Arid . Il989h . 
Disre garding the possibility of viscoelastic effects (jRevl . 
2006), the stress must be a symmetric tensor linear in 
the velocity gradients: 



-,aP 



Ag 



aP 



(148) 
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with A linearly proportional to S a a and B independent 
of S. Since S a a = for an incompressible membrane, 
only B a/3 ^ s enters T^f ■ The symmetries of the stress 
and rate of strain tensors imply that B a ^ lS is symmetric 
under exchange of a and (3, and also 7 and S. In principle, 
we could build B a ^ 1& out of combinations of both of the 
tensors K af3 and g a/3 . However, we will work to leading 
order in an expansion in curvature, in which case 



a(3 -f5 



(g ap g 



(149) 



The constant ( is the dilational viscosity, which does not 
enter the stress for an incompressible membrane since 
S% = 0. Thus, T"£ = 2nS a(3 , and the total membrane 
viscous stress is 



-pa 

vis 



[V a V p + V f3 V a - 2K af3 V n ) t 



0- 



(150) 



The force per unit area acting on the membrane due to 
the internal viscous stress is f V j S = V Q FfL, or 



f . 



2 (Hg af} - 
2^[K a pV c 



KV 13 - 2V n V l3 H 

K a/3 ) V a V n ]t p 

V - V n {AH 2 - 2K)]n, 



(151) 



where the Gauss-Codazzi equations (I104p have been 
used to simplify the final form. Various limit s and 
forms of fyjs have appeared in the literature (| Arid . 1 19891 : 



Arroyo and DeSimond. 120091: 15 orries and Folti n. 1996; 
Fuiitanil . Il997t IScrivenl . Il960l : IWaxmanl . Il984t lYouheil . 



19941 ). Note that Eq. (|15ip leads to a qualitative change 



in the equations for fluid flow even in the simplest case 
of flow on a fixed sphere, for w hich f v j s = jj,CV 2 V ^ + 
KV^)tp, where K is a constant (|Henle et all l2008t h 

It is natural to ask if there are any other terms of 
the same order in gradients of velocity that should be 
included in the surface viscous stress. For example, 
the term F" = /j,2nS7 a V n with "transverse" viscosity 
fi2 has often been used in the theory of capillar y waves 
on surfactant-laden interfaces ([Goodrich! . 1 1 9 6 lh ■ How- 
ev er, this t erm h as recently been shown to be unphysical 
bv lBuzzal (|2002h . Buzza's argument is that F" does not 
obey frame invariance. Frame invariance is a fundamen- 
tal assumption of continuum mechanics that states that 
the relation between stress and strain for a small ele- 
ment of material cannot depend on the motion of the 
element ( Larsonl [l988i : lOldrovdl . Il950h . since at the scale 
of the microstructure, inertial effects are small for typi- 
cal shear rates. In particular, the relation between stress 
and strain for a small element must remain the same 
if the element is subject to a time-dependent rotation. 
Removing the spurious term F„ eliminates unphysical 
results such as negative measured dilatational viscosi- 
ties that have plagued the interpretation of quasi-elastic 
light scattering experiments f or many years ( Buzzal [20021 : 
ICicuta and Hopkinsonl 120041) . 

To see why F" is unphysical, one can show that in a 
rotating frame, the force per unit volume depends on th e 
rate of rotation for a rigid-body rotation ( Buzzal [2002h . 



Perhaps a more direct argument that something is amiss 
with F" is that it predicts a positive di ssipation rate for 
rigid-body rotation of a flat membrane ( Buzzal . [20021 ). If 
we include the [i\ term, then the dissipation rate 



W = / V Q V • F« is dS, 



(152) 



has a term in the integrand proportional to 
M2(V a Vn)(V a V n ), which is nonzero for a rotation 
about an axis lying in the plane of the membrane. Thus 
we conclude that /ii = 0. 

The membrane is also subject to tractions from the 
viscous solvent. The stress tensor for a viscous fluid is 



-pSij + r)(ViUj + VjUi). 



(153) 



For membranes and interfaces, there is no systematic ex- 
pansion that leads to simple approximations analogous to 
slender-body theory, although an approximation similar 
to resistiv e-force theory is sometime s used for illustrative 
purposes ( Cai and Lubenskvl . Il995f h Balancing all the 
forces on a patch of membranes, we find 



V Q (F« ot ) J= 0, 



(154) 



where <7-„- is the three-dimensional viscous stress eval- 
uated at the membrane surface with n the outward- 
pointing normal, and F" ot is the elastic and viscous 
stresses of Eqs. (p^Hj) . p^|) and p5U|) . 

We close this subsection by writing out the equations 
of motion for a membrane which is almost flat. We 
work to lowest order in deflection, and study the evo- 
lution of a ripple with wavelength 2Tr/q. The mem- 
brane shape is give by r(x, y, t) = (x, y, z(x, y, i)), where 
z(x,y,t) = h{t) exp(igx), and complex notation is used 
since we work to linear order in h. The large value of 
the membrane bulk modulus implies that tangential flow 
in the membrane is small, which is why there is no mo- 
tion in the sy-plane in the expression for r(x,y,t). The 
equations of motion for the membrane simplify greatly 
since the membrane flow vanishes, V — 0, and many of 
the geometrical quantities such as K are second order in 
h. The mean curvature is H w V 2 z/2, where to leading 
order V 2 w d 2 /dx 2 + d 2 /dy 2 . The normal component of 
force balance on the membrane is thus 

0+ -a z -) z=Q - K V 4 z = 0. (155) 

To find the stresses af , solve the Stokes equations (|57l - 
I58[) subject to the no-slip boundary conditions u x = 
and u z = hexp(iqx): 

u± = -qzhe iqx T qz (156) 
uf = h(l±qz)e iqx ^ qz (157) 



P 



±2-qqh. 



(158) 



Thus, [af z — <J zz ] z= q = —Arjqhexjpiqx, which together 
with Eq. (|155[) implies 



h = -^h 
4i] 



(159) 
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The ripp le relaxes with the ben ding relaxation 
time (|18l) ( Brochard and Lennonl . Il975[ ) . Note that even 
if 9^sd 3> 1, the membrane viscosity does not enter since 
there are no tangential flows due to the high bulk mod- 
ulus. However, curvature can lead to tangential flows 
for incompressible membranes, as we show in the next 
section. 



E. Illustrative example: buckling and pearling instabilities 
of tubular polymersomes 

In this section we discuss two instabilities of a tubular 
vesicle: pearling and buckling. Pearling arises when the 
membrane is subject to tension, and buckling arises when 
the membrane is under compression. Both instabilities 
consist of a cylinder deforming into a series of bulges, 
but the physical mechanisms are different. Pearling can 
be induc ed by applying a laser tweezer to a tubular l i- 
posome (|Bar-Ziv and MosesL 11994 iBar-Ziv et aZlll998t ). 
The laser tweezer induces a tension, and if the tension 
is sufficiently high, the cylinder can lower its energy by 
deforming into series of bulges. This instability is simi- 
lar to the Rayleigh-Plateau instabili t y of c ylindrical in- 
terfaces (iRavleighl 1 18921: iTomotikal Il935f). and it has 
been extensively studied in liposomes dGoldstein et all . 
Il996bl iGranek and Olamil . 1199,4 iNelson et all . Il995h . 
Buckling of a membrane is ana logous to the buck- 
ling of a rod under tension (jLandau and Lifshitzl . 
Il986t ). Buckling of a sphe r ical m em brane is well- 



p = 3 



studied (iDeuling and Helfrichl. ll976aHbl: IJerikind. Jl977t 
|Petersonfll988lll989l : IZIiong-can and Helfrichll 19871 ). but 
the buckling of a cylindrical membrane has received less 
attention. 

Consider a cylindrical polymersome of radius a. If a 
is of the order of microns, then isn/a* ^ 1, and sur- 
face membrane viscosity cannot be disregarded. We will 
study the instability of the cylinder as a function of the 
pressure p inside the vesicle, with the pressure outside 
taken to be zero. Increasing the pressure increases the 
tension in the membrane and leads to pearling; decreas- 
ing the pressure compresses the membrane and leads to 
buckling. Suppose that in the initial state that there are 
just enough molecules to enclose a cylindrical volume of 
radius a without stretching the membrane. Thus, <j> = 4>q, 
7 = — p-2d — 0, and there is no stretching stress in the 
membrane, F Q = [Eq. (|132p j. However, because the 
membrane is curved, there is a nonzero bending stress 
Fg ond and a corresponding bending force per unit area 



[Eq. (|129p ]. which must balance the internal pressure: 
pn-l- fbcnd = 0. Since H = —l/(2a), normal force balance 
in the initial equilibrium state yields p = po = —k/ (2a 3 ). 

Now imagine changing the pressure to a value p ^= 
Pq. The cylinder stretches or compresses, depending on 
whether p > or p < 0. Since the stretch modulus 
k is high, ka? / n ^> 1, there is little stretch: </>i/<A) = 
(<f> — <A))/0o *C 1- However, the tension 7 is not small [cf. 




FIG. 11 Dimensionless growth rate Xfj,a 2 / 'n vs. dimension- 
less wavenumber qa for the pearling instability for a tubular 
polymersome, in which membrane surface viscosity dominates 
the dissipation. The solid lines have iso/a = 1000 and the 
dashed line has iso/a = 00. 



discussion after Eq. (|132[l ]. and 

7 =p + n/(2a 3 ). 



(160) 



Our task is to determine if sufficiently large (positive or 
negative) tension in the membrane destabilizes the cylin- 
drical shape. 

To proceed, parameterize the vesicle in cylindrical co- 
ordinates, r((f, z, t) = p(a + h(z 7 1)) + (z + w(z 7 t))i. Note 
that the deformed shape is assumed axisymmctric. To 
first order in the displacements h and w, the mean cur- 
vature is H ps — (l/a+ h/a? + h")/2, and the Gaussian 
curvature is K sa —h"/a. The evolution of the shape of 
the membrane is governed by normal force balance, 



0. 



(161) 



njcr^(p = a,z,t)hj is the normal- 



n • (f vis + f b 

where <r± n \ p = a 
normal component of the stress tensor for the solvent 
outside (+) and inside (— ) the vesicle evaluated at the 
membrane, to leading order in h. 

We suppose the deformation is sinusoidal, h(z,t) cx 
cxpi(oz). Again, since ka? f k, 3> 1, the change in density 
is small and we may use the condition of continuity for an 
incompressible membrane, Eq. (|147p . to relate the shape 
to the flow V z = w in the membrane. To linear order, 



igl^ 2 



-h/a. 



(162) 



The internal forces per unit area acting on the membrane 
arc thus 



n 



(fv 



Ibcnd 



f s ) re -2^h/a 2 - £{q)h, (163) 



where £{q) — k(o 4 o 2 



pq 2 a 2 



1 — p)/a , and we have 



introduced the dimensionless pressure p = pa? j k. 

The solvent forces acting on the membrane are deter- 
mined by solving Stokes equations (J57] [58]) for the flow 



inside and outside the tube, and then calculating 
the stress a^ n . Note that the imposed pressure p en- 
ters the stress tensor through the dynamic pressure P: 
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P ± Sij + r)(diU^ + djuf), where P = p when 

h = 0, and P + = when h = 0. Since the Stokes equa- 
tions are linear, the stress at the membrane has the form 
( a nn ~ <J nn)p=a = —rjK^ 1 h + p, where the "dynamical 
factor" A depends only on the wavenumber q and radius 
a. The Stokes equations in an axisymmetric geometry 
are readily solved u sing the stream func t ion where 
(ip^z) ( Happel and Brenner . 1965b ). Here 



u 



we on ly quote the result and refer to ( Goldstein et all 
199Gb) for the details. With the no-slip boundary condi- 
tions at the membrane, u^(a) — h and uf(a) = V z , the 
dynamical factor is 



A 



a [qa(I 2 - If) - 2I h] [ga(K 2 - K 2 ) + 2KqK^ 

2 1 X K X + I K + qa(hK - I K{) 

(164) 

where /„ and K n are modified Bessel functions of order 
n evaluated at qa. Note that this expression correct s 
an error in equation (D.9) of ([Goldstein et all Il996blh 
which gave a dynamical factor which is too small. 

Balancing internal membrane forces and the solvent 
force leads to h = Xh, where 



k q 4 a A + pq 2 a 2 + 1 



P 



Ha? 2+(a/A)(a/e ST) ) 



(165) 



The 



sign 



of A 



is determined by the sign of £{q) = 
1 -p)/a 4 . When p > n/a 3 , £ > 
and the cylinder is unstable. By Eq. p60[) , this pressure 
corresponds to a tension 7 = 3 K/(2a 3 ), the critical ten 



K(q A a A + pq 2 a 2 



sion for the pearling instability ( Goldstein et a/.l . ll996bh . 
Figure [Tl] shows the growth rate A for the pearling in- 
stability for various pressures and Isd/cl = 1000. Also 
shown for comparison is the growth rate for p = 3«/a 3 
and ^sd/o = 00. Note how A, which arises from the dy- 
namics of the incompressible solvent, forces the growth 
rate to vanish at q = 0, even when £sd/i 3* 1- 

When p < -2(1 + V2)n/a 3 sa -4.828K/a 3 , the cylin- 
der is also unstable. Since negative pressure corresponds 
to compressive stress, this instability is analogous to the 
Euler buckling instability of a rod. Figure Q2] shows the 
growth rate for the buckling instability. Note the con- 
trast with the pearling instability. In the buckling insta- 
bility, when the pressure is just sufficient to induce the 
instability, only a narrow band of wavenumber centered 
about a nonzero wavenumber [q*a = (1 + V^)] has pos- 
itive growth rate. In the pearling instability, the band 
extends from q = to a finite value, and therefore there 
is no limit to the wavelength of the unstable modes. 




FIG. 12 Dimensionless growth rate Xfj,a 2 / ' n vs. dimension- 
less wavenumber qa for the buckling instability for a tubular 
polymersome, in which membrane surface viscosity dominates 
the dissipation. The solid lines have iso/a = 1000 and the 
dashed line has £so/a = 00. 



such as curvature. We have also seen how naturally these 
geometric ideas fit with the variational point of view. In 
the case of filaments, the variation of the energy led di- 
rectly to the constitutive relation between moments and 
the rate of rotation of the directors, or strain. In the 
case of membranes, variation using auxiliary variables 
provided a quick route to the elastic forces per area, with 
the added bonus of the explicit form of the stress tensor 
of the membrane. Geometric ideas are also crucial for un- 
derstanding viscous flow on a deforming curved surface. 
This study of polymers and membranes remains an ac- 
tive area of research, and the purely mechanical problems 
considered here provide a solid foundation for probing the 
dynamics of fluctuating polymers and membranes. 
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APPENDIX A: Coordinate transformations 



V. OUTLOOK 

The aim of this review has been to give a practical 
overview of the geometrical tools necessary to formu- 
late the equations of motion of polymers and membranes. 
Equilibrium and dynamic properties of thin filaments and 
membranes are best described using geometrical ideas 



In this appendix we describe how the components of 
tensors transform under a change of coordinates. In gen- 
eral, a tensor is a quantity that transforms linearly under 
a change of coordinates. We begin with the metric tensor 
as a concrete example. 

Consider two systems of coordinates related by = 
f a (^ 1 ,^ 2 ). From the chain rule, the metric tensor g a p in 
the new coordinate system is related to metric in the old 
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coordinates by a linear transformation law: 
dr dr d£ 5 



(Al) 



The differentials have the opposite transformation law, 

(A2) 

which ensures that 

gapd^df = g a pdCd^. (A3) 

Quantities that transform linearly under a change of co- 
ordinates are tensors. The tensor transformation law for 
g a /3 reflects the fact that the distance ds 2 is a geometric 
quantity, independent of the choice we make for coordi- 
nates. 

It is traditional to use upper or lower indices to indicate 
what kind of transformation law a tensor has. Thus, a 
tensor T^^'" transforms homogeneously under a change 
of coordinates as 



rpOLiO.2... 



n 



<9£ 71 <9£ 72 d£fr d£fr 5l5: 



(A4) 



Tensors with upper indices are "contravariant," and ten- 
sors with lower indices are "covariant." 

As mentioned in PV[ we use raised indices to denote 
the inverse of the metric tensor: 



g^g lP = S°p. 



(A5) 



Viewing Eq. (|A1|) as a matrix equation and inverting 
yields 



9 



a/3 



dp d£ s 9 ' 



(A6) 



that is, g a @ transforms as a contravariant tensor. By 
using the metric tensor, we can create a covariant tensor 
from a contravariant one, and vice-versa. For example, 
W a = g a pW and W a = g al3 W . Geometrically, we 
can think of W a as the components of the vector field 
W = W a t a with respect to a dual basis t Q defined by 



I 



0- 



(A7) 



In other words, W = W a t a = W a t a . The dual basis 
vectors t Q are like reciprocal lattice vectors, and also obey 
the rule about raising and lowering indices: 



9 



a/3. 



(A8) 



Recall that the covariant derivative was defined in fllVI 
to be a geometric object, independent of coordinates once 
the direction t a is chosen. The quantity (V Q V r )' 3 = f 3 • 
9 Q V is a tensor since 9 Q V and t^ 3 are both tensors. Note 
however that the Christoffel symbol V^ g is not a tensor. 



The importance of these linear transformation laws 
is that they allow us to easily construct coordinate- 
invariant expressions. A quantity is coordinate invariant 
when there is an equal number of upper and lower indices 
that are summed over. For example, the mean curvature 
H = K a a and the divergence W a V a are both invariant 
under coordinate changes. Likewise, the Gaussian curva- 
ture det K a g is invariant under coordinate changes, since 
the upper and lower indices have associated with them 
transformation matrices that are inverses. On the other 
hand, g — detg a p is not coordinate-invariant since it has 
two lower indices and no upper indices, and therefore gets 
two factors of the Jacobian matrix \d(^, ^ 2 )/9(^" 1 , £ 2 )| 
under a coordinate change. These factors are precisely 
what is needed to make dS = */gd£ d£ 2 invariant under 
coordinate change. 



APPENDIX B: Green's theorem 

Consider Stokes theorem for a tangent vector field A 
A a t a on a region D of a surface, 



n ■ V x Ad,? 



A ■ ds, 



(Bl) 



i)D 



where ds is the infinitesimal tangent vector along bound- 
ary, and V denotes the three-dimensional gradient in 
Cartesian coordinates. Note that V is the covariant 
derivative for Euclidean space in standard coordinates. 
Recall that the normal component of the curl has deriva- 
tives along the tangent direction only: 

n • V x A = (nxV)-A 
~ti x t 2 



V9 



1 

71 
i 



x V • A 
[ta(ti-V)-ti(t 2 .V)]-i4°t a 



(ViA 2 - V 2 Ai) 



= — (diA 2 
V9 



Since A • ds = A a d^ a , we have shown that 
[ e^VoApdS = j>A a dC, 



(B2) 



(B3) 



,-21 



where we have introduce the two-dimensional antisym- 
metric tensor with components e 12 = — 
e ii = e 22 _ o. Defining e Q/3 = g aj gf3S£ lS 
and A a = ep a B°, we also have 



1/y/g and 
(note ei2 = y/g) 



V a B a dS 



je^B^dC- 



(B4) 



A short calculation shows that e a /3B a d^ — B • pds, 
where B = B a t a ; p is the outward-pointing unit vec- 
tor normal to the curve dD defining the region of in- 
tegration and lying in the tangent plane, p = t Q d£ a x 
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n/Vdp-yd^d^ = dr/ds x n; and ds = y/g^d^d^ is 
measure of arclength along the curve. Therefore we may 
write 

/ V a B a dS = I B a p a ds (B5) 

J D JdD 

where p a ds — e Q ^d^. 
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